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ABSTRACT 

We  present  a  study  of  the  [C  ii]  157.74  /mi  fine-structure  line  in  a  sample  of  15  ultraluminous  infrared  (IR) 
galaxies  (IR  luminosity  UIR$ TO12  L0;  ULIRGs)  using  the  Long  Wavelength  Spectrometer  (LWS)  on  the 
Infrared  Space  Observatory  (ISO).  We  confirm  the  observed  order  of  magnitude  deficit  (compared  to  normal 
and  starburst  galaxies)  in  the  strength  of  the  [C  ii]  line  relative  to  the  far-infrared  (FIR)  dust  continuum 
emission  found  in  our  initial  report,  but  here  with  a  sample  that  is  twice  as  large.  This  result  suggests  that  the 
deficit  is  a  general  phenomenon  affecting  4  out  of  5  ULIRGs.  We  present  an  analysis  using  observations  of 
generally  acknowledged  photodissociation  region  (PDR)  tracers  ([C  n],  [O  i]  63  and  145  /on,  and  FIR  contin¬ 
uum  emission),  which  suggests  that  a  high  ultraviolet  flux  Go  incident  on  a  moderate  density  n  PDR  could 
explain  the  deficit.  However,  comparisons  with  other  ULIRG  observations,  including  CO  (1-0),  [C  i]  (1-0), 
and  6.2  pm  polycyclic  aromatic  hydrocarbon  (PAH)  emission,  suggest  that  high  Go / n  PDRs  alone  cannot 
produce  a  self-consistent  solution  that  is  compatible  with  all  of  the  observations.  We  propose  that  non-PDR 
contributions  to  the  FIR  continuum  can  explain  the  apparent  [C  ii]  deficiency.  Here,  unusually  high  G0  and/ 
or  n  physical  conditions  in  ULIRGs  as  compared  to  those  in  normal  and  starburst  galaxies  are  not  required 
to  explain  the  [C  ii]  deficit.  Dust-bounded  photoionization  regions,  which  generate  much  of  the  FIR  emission 
but  do  not  contribute  significant  [C  ii]  emission,  offer  one  possible  physical  origin  for  this  additional  non- 
PDR  component.  Such  environments  may  also  contribute  to  the  observed  suppression  of  FIR  fine-structure 
emission  from  ionized  gas  and  PAHs,  as  well  as  the  warmer  FIR  colors  found  in  ULIRGs.  The  implications 
for  observations  at  higher  redshifts  are  also  revisited. 

Subject  headings:  galaxies:  high-redshift  —  galaxies:  ISM  —  galaxies:  starburst  —  infrared:  galaxies  — 

ISM:  lines  and  bands 


1.  INTRODUCTION 

In  a  previous  study  (Luhman  et  al.  1998,  hereafter  Paper 
I),  we  reported  measurements  of  the  157.74  pm  2Pi/2~2P\/2 
fine-structure  line  of  C+  in  a  sample  of  seven  ultraluminous 
infrared  (IR)  galaxies  (IR  luminosity  L\\\  >  1012  LQ; 
ULIRGs)  using  the  Long  Wavelength  Spectrometer  (LWS) 
on  the  Infrared  Space  Observatory  (ISO).  The  observations 
presented  in  Paper  I  identified  for  the  first  time  a  [C  ii]  deficit 
in  ULIRGs  relative  to  the  far-infrared  (FIR)  dust  contin¬ 
uum,  whereby  the  [C  ii]  flux  is  about  an  order  of  magnitude 
lower  than  that  seen  from  nearby  normal  and  starburst  gal¬ 
axies.  As  an  important  coolant  of  atomic  and  molecular 
interstellar  gas  exposed  to  ultraviolet  (UV)  photons,  com¬ 
monly  referred  to  as  photodissociation  regions  or  PDRs 
(Tielens  &  Hollenbach  1985),  the  [C  ii]  158  pm  line  is  often 
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the  single  brightest  line  in  the  spectrum  of  most  galaxies, 
providing  0. 1%— 1%  of  the  total  FIR  luminosity  (see,  e.g., 
Stacey  et  al.  1991  and  references  therein).  However,  based 
on  the  Paper  I  observations,  the  [C  ii]  contribution  in 
ULIRGs  is  more  typically  0.01%  0. 1%  of  the  total  FIR 
luminosity. 

The  unexpected  low  [C  ii]  luminosity  result  could  have 
important  consequences  for  future  spectroscopic  surveys  of 
protogalaxies  from  the  ground  or  space  (e.g.,  Stark  1997; 
Blain  et  al.  2000).  The  ULIRGs  in  this  study  are  local  exam¬ 
ples  of  merging  galaxies  (see  Sanders  &  Mirabel  1996  for  a 
review).  Since  protogalaxies  at  high  redshift  (z  =  1-10)  are 
believed  to  form  from  the  merging  of  smaller  structures 
(e.g.,  Mobasher,  Rowan-Robinson,  &  Georgakakis  1996; 
van  den  Bergh  et  al.  1996),  they  could  resemble  the  nearby 
ULIRGs.  As  the  brightest  line  in  the  spectrum  of  most  gal¬ 
axies  and  particularly  bright  in  low-metallicity  galaxies 
(Madden  2000),  the  [C  ii]  line  is  a  potentially  powerful 
spectroscopic  tracer  of  protogalaxies,  where  for  redshifts 
z  =  2-6,  the  [C  ii]  158  pm  line  shifts  to  submillimeter  wave¬ 
lengths.  Predictions  of  the  detection  rate  of  the  [C  ii]  line  as 
a  signpost  of  high-redshift  galaxies  must,  however,  take  into 
account  the  [C  n]  deficit  found  in  Paper  I. 

In  Paper  I,  we  did  not  have  sufficient  data  in  hand  to 
explore  in  depth  the  source  of  the  [C  ii]  deficiency  in 
ULIRGs.  However,  since  the  publication  of  Paper  I,  we 
have  obtained  additional  FIR  spectra — increasing  our  sam¬ 
ple  size  to  15  ULIRGs — to  confirm  the  [C  ii]  deficit  and 
readdress  its  origin.  In  some  cases,  we  combined  our  added 
spectra  with  archival  [C  ii]  data  from  other  ISO  programs. 
In  addition  to  the  added  FIR  spectra,  several  new  studies 
have  been  performed  that  provide  useful  clues  regarding  the 
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nature  of  the  [C  n]  emission  in  ULIRGs.  Satyapal  et  al. 
(2003)  and  S.  Satyapal  et  al.  (in  preparation)  performed 
detailed  analyses  of  the  ISO  FIR  spectra  of  the  luminous 
galaxy  Arp  299  and  applied  these  analyses  to  a  discussion  of 
the  FIR  line  emission  in  ULIRGs  in  general;  Fischer  et  al. 
(1999,  2001)  reported  on  the  emission  and  absorption  line 
characteristics  of  the  FIR  spectra  of  a  small  sample  of  IR- 
bright  galaxies,  including  the  ULIRG  Arp  220,  and 
modeled  the  photoionization  of  dust-bounded  nebulae;  and 
Genzel  et  al.  (1998),  Lutz  et  al.  (1998),  Rigopoulou  et  al. 
(1999),  Dudley  (1999),  Tran  et  al.  (2001),  and  C.  C.  Dudley 
et  al.  (in  preparation)  studied  the  polycyclic  aromatic 
hydrocarbon  (PALI)  emission  in  luminous  IR  galaxies. 
Several  of  these  studies  tend  to  suggest  the  probable  impor¬ 
tance  of  an  additional  component  of  dust  emission  not 
traced  by  PDR  line  diagnostics. 

In  related  studies,  Malhotra  et  al.  (2001)  analyzed  the 
FIR  spectra  of  a  sample  of  60  normal,  star-forming  galaxies. 
For  the  same  sample  of  galaxies,  Helou  et  al.  (2001)  com¬ 
pared  the  [C  ii]  line  emission  to  the  5-10  /mi  band  aromatic 
feature  emission  and  found  evidence  for  a  decrease  in  the 
importance  of  PAHs  and  small  grains  relative  to  large 
grains  with  increasing  radiation  field  both  in  gas  heating 
and  in  dust  cooling.  For  those  FIR-bright  galaxies  in  the 
Malhotra  et  al.  sample  that  display  a  [C  ii]  deficit,  Malhotra 
et  al.  attributed  the  decrease  in  the  [C  n]  to  FIR  flux  to 
PDRs  with  a  high  ratio  of  incident  UV  radiation  G0  to  gas 
density  n  based  on  the  values  of  G0  and  n  derived  from 
plane-parallel  PDR  models.  In  a  similar  study,  however, 
Negishi  et  al  (2001)  found  that  the  ratio  of  6),  to  n  remains 
constant  in  their  sample  of  34  galaxies  and  surmised  instead 
that  galaxies  with  low  ratios  of  [C  ii]  to  FIR  flux  could  be 
interpreted  in  terms  of  a  decrease  in  the  ionized  component 
of  the  [C  ii]  emission  relative  to  the  FIR  intensity.  Using  the 
CO  (1-0)  and  [C  i]  (1-0)  line  emission,  several  authors 
(Genzel  &  Cesarsky  2000;  Gerin  &  Phillips  2000)  have  also 
recently  derived  PDR  physical  conditions  for  ULIRGs.  The 
derived  physical  parameters  for  these  various  PDR  studies 
vary  significantly  and  show  that,  depending  on  which  lines 
are  used,  PDR  analyses  can  give  different  results. 

In  this  paper,  we  readdress  the  [C  ii]  deficit  in  ULIRGs, 
combining  the  results  of  these  recent  ULIRG  studies  with 
our  own  ISO  LWS  observations  of  [C  ii]  as  well  as  archival 
data  from  other  ISO  programs  for  an  enlarged  sample  of 
ULIRGs.  Along  with  the  [C  ii]  data,  we  also  report  [O  i] 
63  /an  measurements  and  [O  i]  145  /mi  upper  limits  in  a  subset 
of  the  sample  ULIRGs.  Like  the  [C  ii]  158  /mi  line,  the  [O  i] 
63  //m  line  and,  to  a  much  lesser  degree,  the  [O  i]  145  /mi  line 
are  strong  cooling  lines  in  Galactic  and  extragalactic  PDRs 
(e.g.,  Tielens  &  Hollenbach  1985;  Wolfire,  Tielens,  & 
Hollenbach  1990).  The  combined  [C  ii]  and  [O  i]  data  set, 
along  with  the  FIR  emission,  allows  us  to  derive  estimates  of 
Go  and  n  using  the  predictions  of  PDR  models  and  assuming 
that  the  [C  ii],  [O  i],  and  FIR  emission  each  emanate  from  the 
same  effective  PDR  surface  area.  In  order  to  treat  the  PDR 
emission  in  a  self-consistent  way,  we  compare  the  [C  ii],  [O  i], 
and  FIR  emission  to  other  PDR  diagnostics,  namely  the  CO 
(1-0),  C  i  (1-0),  and  6.2  /mi  PAH  flux,  and  detennine  the 
proper  application  of  one-dimensional,  plane-parallel  PDR 
models  to  ULIRG  observations.  Our  analysis  suggests  the 
presence  of  a  contributing  component  of  non-PDR  FIR 
emission  to  the  total  FIR  luminosity  in  ULIRGs.  We  discuss 
the  consequences  of  this  result  on  the  usefulness  of  [C  ii]  as  a 
tracer  of  high-redshift  protogalaxies. 


2.  OBSERVATIONS  AND  RESULTS 

In  Table  1,  we  summarize  the  ULIRG  observations.  The 
ISO  programs  from  which  we  extracted  [C  n]  data  include 
the  ISO  LWS  consortium  “Central  Programme”  on 
IR-bright  galaxies  (IRBGALS  and  IRBGALS2),  which 
provided  the  data  for  Paper  I;  the  ISO  OH  megamasers  pro¬ 
gram  (OHMEGA);  our  follow-up  ISO  discretionary-time 
program  on  the  [C  ii]  emission  in  ULIRGs  (CII_ULIRG); 
and  the  Max-Planck-Institute  fur  Extraterrestrische  Physik 
extragalactic  program  (MPEXGAL2).  We  also  extracted 
archival  data  from  several  other  ISO  programs  (LWS_PVL 
IRB_HD,  and  SF_GLX_A),  which  were  found  to  contain 
[C  ii]  spectra.  In  general,  the  IRBGALS  and  CII_ULIRG 
programs  contain  observations  of  the  [C  n]  158  /an  line 
only,  whereas  the  MPEXGAL2  program  contains  both 
[C  ii]  158  /an  and  [O  i]  63  /an  observations.  These  line  obser¬ 
vations  used  the  ISO  LWS  medium-resolution  (AA  =  0.6 
/an  at  158  /mi  or  R  ~  200)  grating  mode  (Clegg  et  al.  1996). 
In  the  cases  of  Mrk  231,  Mrk  273,  Arp  220,  and  IRAS 
17208—0014,  we  obtained  full  43-197  /an  spectra,  from 
which  we  extracted  flux  measurements  of  the  [C  ii]  158  /an 
and  [O  i]  63  and  145  /an  lines. 

Pipeline  processing  (pipeline  10)  of  the  raw  ISO  LWS 
spectra  calibrated  the  wavelength  and  flux,  removed 
cosmic-ray  glitches,  subtracted  dark  current,  and  corrected 
for  instrumental  responsivity  variations  (Swinyard  et  al. 
1996).  We  performed  subsequent  data  processing,  including 
co-addition,  line  profile  extraction,  and  line  flux  and  uncer¬ 
tainty  measurement,  using  the  ISO  Spectral  Analysis 
Package  (ISAP9;  Sturm  et  al.  1998).  For  sources  for  which 
multiple  observations  were  taken  from  different  programs, 
the  line  fluxes  are  the  weighted  means  with  propagated 
errors.  Table  2  contains  the  enlarged  sample  of  ULIRGs 
and  their  corresponding  [C  ii]  line  fluxes  or  upper  limits. 
The  [C  ii]  fluxes  here  using  pipeline  10  processing  compare 
well  with  the  pipeline  5  fluxes  reported  in  Paper  I,  although 
because  of  the  co-addition  and  use  of  pipeline  10,  the  errors 
have  decreased  on  average.  The  present  study  increases  the 
number  of  ULIRGs  in  our  sample  from  7  to  15.  (We  note, 
however,  that  LIR  is  slightly  less  than  1012  Le  for  two  gal¬ 
axies  in  our  sample,  IRAS  19254—7245  and  IRAS 
23128—5919;  these  galaxies  are  not  strictly  defined  as 
ULIRGs.  To  within  the  uncertainties  of  the  Ljr  estimates, 
however,  these  two  galaxies  are  indistinguishable  from  1012 
Le  galaxies,  so  we  include  them  in  our  ULIRG  sample.)  In 
addition.  Table  2  lists  the  [O  i]  63  and  145  /an  line  flux  meas¬ 
urements  and  upper  limits  where  available.  The  absolute 
flux  calibration  error  is  ±25%  and  was  determined  by 
comparing  overlapping  data  in  neighboring  detectors 
from  the  full-grating  scans.  All  sources  in  this  study  were 
spatially  unresolved,  with  the  ISO  beam  (69"  at  158  /an) 
encompassing  all  of  the  IR  emission. 

As  shown  in  Table  2,  the  [C  ii]  line  has  been  detected  in  12 
of  the  1 5  ULIRGs  in  our  sample.  Of  the  nine  ULIRGs  for 
which  we  have  63  /an  spectra,  the  [O  i]  63  /an  line  is  seen  in 
absorption  in  Arp  220  (Fischer  et  al.  1999),  and  in  emission 
in  Mrk  231,  Mrk  273,  IRAS  17208-0014,  IRAS 
20551-4250,  and  IRAS  23365±3604.  We  measured  [O  i] 
145  /an  upper  limits  for  Mrk  231,  Mrk  273,  Arp  220,  and 
IRAS  17208-0014. 

In  Figure  1,  we  plot  the  [C  ii]  line  flux  against  the  FIR  flux 
for  the  ULIRGs  in  our  sample  and  for  an  assortment  of  nor¬ 
mal  and  starburst  galaxies  from  previous  studies  (Paper  I; 


TABLE  1 


Observing  Log 


Source 

(1) 

«2000 

(2) 

62000 

(3) 

750  Program 
(4) 

AOT 

(5) 

Observation  Date 
(6) 

IRAS  05189-2524 . 

0521  01.4 

-25  21  44.9 

CIIJJLIRG 

LWS02 

1998  Apr  3 

IRAS  12071-0444 . 

12  09  45.3 

-05  01  13.3 

IRBGALS 

LWS02 

1996 Jul 19 

Mrk231 . 

12  56  14.2 

+56  52  24.9 

LWS_PV1 

LWS01 

1996  Jan  7 

12  56  14.2 

+56  52  24.9 

IRBGALS 

LWS01 

1996  May  15 

12  56  14.2 

+56  52  24.9 

IRBGALS 

LWS01 

1997 Jul 11 

Mrk  273 . 

13  4442.1 

+55  53  13.2 

MPEXGAL2 

LWS02 

1996  Apr  25 

13  4442.1 

+55  53  13.2 

IRBGALS 

LWS02 

1996  May  19 

134441.6 

+55  53  18.7 

OHMEGA 

LWS01 

1997 Jul 22 

IRAS  15206+3342 . 

152238.1 

+33  31  33.3 

IRBGALS 

LWS02 

1997  Sep  5 

IRAS  15250+3609 . 

15  26  59.3 

+35  58  37.2 

IRBGALS 

LWS02 

1996  Sep  20 

Arp  220 . 

15  34  57.2 

+23  30  11.2 

IRBGALS 

LWS01 

1996  Aug  20 

15  34  57.2 

+23  30  11.2 

IRBGALS 

LWS02 

1997 Jul 18 

15  34  57.2 

+23  30  11.2 

IRB_HD 

LWS01 

1997  Aug  17 

IRAS  17208-0014 . 

172321.9 

-0017  01.0 

MPEXGAL2 

LWS02 

1996  Mar  5 

1723  21.9 

-00  17  01.0 

OHMEGA 

LWS01 

1997  Aug  26 

IRAS  19254-7245 . 

193121.6 

-72  39  24.9 

MPEXGAL2 

LWS02 

1996  Sep  11 

IRAS  19297-0406 . 

193222.1 

-04  00  02.2 

CIIJJLIRG 

LWS02 

1998  Mar  21 

IRAS  20100-4156 . 

20  13  29.8 

-41  47  34.7 

OHMEGA 

LWS02 

1996  Nov  9 

2013  29.8 

-41  47  34.7 

CIIJJLIRG 

LWS02 

1998  Mar  21 

20  13  29.8 

-41  47  34.7 

IRBGALS2 

LWS02 

1998  Mar  23 

IRAS  20551-4250 . 

20  58  26.9 

-42  39  06.2 

MPEXGAL2 

LWS02 

1996  Apr  20 

20  58  26.9 

-42  39  06.2 

CIIJJLIRG 

LWS02 

1997  Nov  9 

IRAS  22491-1808 . 

22  51  49.3 

-17  52  24.1 

IRBGALS 

LWS02 

1996  May  21 

IRAS  23128-5919 . 

23  15  47.0 

-59  03  16.9 

MPEXGAL2 

LWS02 

1996  Apr  20 

23  15  47.0 

-59  03  16.9 

CIIJJLIRG 

LWS02 

1997  Nov  27 

IRAS  23365+3604 . 

233901.2 

+3621  10.0 

SF_GLX_A 

LWS02 

1997  Dec  6 

23  39  01.2 

+3621  10.0 

CIIJJLIRG 

LWS02 

1997  Dec  7 

Note. — Col.  (1):  source  name;  col.  (2):  right  ascension,  epoch  2000,  in  hours,  minutes,  and  seconds;  col. 
(3):  declination,  epoch  2000,  in  degrees,  arcminutes,  and  arcseconds;  col.  (4):  ISO  program  name;  col. 
(5):  AOT  =  Astronomical  Observation  Template;  LWS01  =  medium-resolution  ( R  ~  200)  spectrum  covering  the 
full  43-197  /an  range  of  the  LWS,  LWS02  =  medium-resolution  line  spectrum  (3.4  fim  width  in  the  case  of  the 
[C  n]  line);  col.  (6):  ISO  observation  date. 


TABLE  2 


Sample  Properties  and  Measured  Line  Fluxes 


Source 

(1) 

D 

(Mpc) 

(2) 

log  (LlR/LQ) 
(3) 

7+TR 

(10_13ergs  cm-2  s-1) 
(4) 

-F[Cn]  158 

( 1 0— 13  ergs  cm-2  s-1) 
(5) 

-f[Oi]63 

( 1 0— 13  ergs  cm-2  s-1) 
(6) 

*ip.]  145 

(10-13  ergs  cm'2  s_1) 

(7) 

TjCn]  158/7+TR 
(8) 

IRAS  05189-2524 . 

170 

12.10 

5870 

1.33  +  0.16 

2.27  x  10“4 

IRAS  12071-0444 . 

514 

12.29 

1110 

<0.6 

<5.41  x  10-4 

Mrk  231 . 

169 

12.54 

14200 

3.51  +  0.14 

2.89  +  0.53 

<0.42 

2.47  x  10~4 

Mrk  273  . 

151 

12.11 

9760 

5.24  +  0.16 

4.70  +  0.67 

<0.30 

5.37  x  10-4 

IRAS  15206+3342 . 

498 

12.18 

812 

1.87  +  0.17 

2.30  x  10  3 

IRAS  15250+3609 . 

222 

12.00 

3110 

<1.6 

<5.14  x  10“4 

Arp  220 . 

73 

12.19 

47900 

9.99  +  0.24 

—6.35  ±  0.73“ 

<0.94 

2.09  x  10-4 

IRAS  17208-0014 . 

171 

12.32 

14500 

6.86  +  0.18 

6.53  +  0.96 

<0.60 

4.73  x  10~4 

IRAS  19254-7245 . 

247 

11.99 

2550 

2.71  +  0.33 

<2.5 

1.06  x  10  3 

IRAS  19297-0406 . 

343 

12.36 

3270 

2.23  +  0.11 

6.82  x  10“4 

IRAS  20100-4156 . 

518 

12.56 

2350 

0.65  ±  0.09 

<4.0 

2.77  x  10~4 

IRAS  20551-4250 . 

171 

12.00 

5410 

3.51  +  0.13 

5+1 

6.54  x  10-4 

IRAS  22491-1808 . 

311 

12.10 

2320 

<1.1 

<4.74  x  10“4 

IRAS  23128-5919 . 

178 

11.95 

4870 

5.58  +  0.26 

<5.0 

1.15  x  10-3 

IRAS  23365+3604 . 

258 

12.13 

3360 

1.32  +  0.10 

1.26  +  0.39 

3.93  x  IQ'4 

Note. — Col.  (1):  source  name;  col.  (2):  distance  for  Ho  =  75  km  s-1  Mpc-1  and  qo  =  0.5;  col.  (3):  8-1000  /an  luminosity  as  measured  by  IRAS  using  the 
prescription  given  in  Table  1  of  Sanders  &  Mirabel  (1996)  [L(8— 1000  //m)  =  5.6  x  105£)JIpc(13.48S'i2  +  5. 16S25  +  2.58560  +  Sroo),  where  5  is  the  flux  density 
in  Jy];  col.  (4):  42.5-122.5  /an  emission  as  measured  by  IRAS  Faint  Source  Catalog  using  the  prescription  given  in  Table  1  of  Sanders  &  Mirabel  (1996) 
[+fir  =  1.26  x  10_11(2. 58560  +  5mo)  ergs  cm-2  s-1]  or  as  measured  by  750  LWS  using  ISAP  (see  text  for  details);  cols.  (5)-(7):  [C 11]  1 58  fim,  [0 1]  63  /an,  and 
[O  1]  145  fim  line  fluxes  (flux  upper  limits  were  calculated  assuming  a  Gaussian  line  with  the  effective  instrumental  profile  and  a  3  cr  amplitude);  col.  (8):  flux 
ratio  of  the  [C 11]  1 58  fim  line  to  FIR  continuum  emission. 

“  The  line  was  seen  in  absorption,  as  denoted  by  the  minus  sign. 
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Fig.  1. — [C  ii]  158  //m  line  flux  observed  with  the  LWS  vs.  the  FIR  flux 
for  15  ULIRGs  compared  with  a  sample  of  normal  and  starburst  galaxies 
(Paper  I;  Colbert  et  al.  1999;  Stacey  et  al.  1999;  Bradford  et  al.  1999; 
J.  Fischer  et  al.,  in  preparation).  In  the  symbol  key,  the  ULIRGs  are  listed 
in  order  of  luminosity  from  top  to  bottom  and  are  denoted  by  filled 
symbols,  while  the  lower  luminosity  galaxies  are  shown  by  open  symbols. 
The  dashed  lines  mark  the  regime  typical  of  normal  and  starburst  galaxies. 


Bradford  et  al.  1999;  Colbert  et  al.  1999;  Stacey  et  al.  1999; 
C.  C.  Dudley  et  al.,  in  preparation;  J.  Fischer  et  al.,  in  prepa¬ 
ration).  Galaxies  are  identified  by  name  according  to  the 
key  in  the  figure.  The  comparison  sample  of  normal  and 
starburst  galaxies  is  limited  to  sources  with  ISO  6.2  /im  spec¬ 
troscopy,  where  the  beam  includes  the  whole  galaxy  or  the 
same  region  as  the  LWS  beam.  Among  the  comparison  sam¬ 
ple  of  galaxies,  four  galaxies — NGC  6240,  NGC  4945, 
NGC  1068,  and  Circinus — with  optical,  X-ray,  or  gamma 
ray  signatures  of  active  galactic  nuclei  (AGNs)  are  included 
(e.g.,  Done,  Madejski,  &  Smith  1996;  Vignati  et  al.  1999). 

We  find  that  the  enlarged  sample  of  ULIRGs  presented 
here  confirms  and  reinforces  the  result  presented  in  Paper  I; 
that  is,  relative  to  the  FIR  continuum,  the  [C  ii]  flux  from 
ULIRGs  is  approximately  an  order  of  magnitude  lower 
than  that  seen  from  nearby  normal  and  starburst  galaxies 
(as  well  as  the  AGNs).  Of  the  15  ULIRGs,  only  three 
objects— IRAS  15206+3342,  IRAS  19254-7245,  and  IRAS 
23128-5919 — display  values  of  /r[Cn]//+iR  within  the  range 
0. 1%— 1%  indicative  of  virtually  all  normal  and  starburst 
galaxies  (Crawford  et  al.  1985;  Stacey  et  al.  1985,  1991; 
Wright  et  al.  1991).  The  FjCn]  to  Ffir  ratio  for  the  remaining 
80%  of  the  sample  ULIRGs  is  less  than  0.1%,  with  a  mean 
value  of  0.041%  excluding  sources  with  upper  limits.  The 
median  of  all  ULIRGs  in  our  sample  is  0.039%,  where  for 
those  ULIRGs  with  [C  n]  upper  limits  we  use  the  2  a  value. 
This  median  value  is  13%  of  the  median  /jCll|  to  Ffir  ratio 
(0.31%)  for  the  60  normal  galaxies  in  the  Malhotra  et  al. 
(2001)  sample,  again  using  the  2  a  value  for  those  normal 
galaxies  with  [C  ii]  upper  limits  only.  (The  ratio  of  the  aver¬ 
age  value  of  our  sample  [0.062%]  and  the  Malhotra  et  al. 
sample  [0.31%]  is  20%.)  For  both  the  Malhotra  et  al.  sample 
and  our  sample,  Arp  220  shows  the  lowest  detected  F|ClI]  to 
CpiR  ratio  of  0.021%.  NGC  4418  has  the  lowest  upper  limit 
of  0.0051%  (Malhotra  et  al.  2001)  of  any  galaxy  in  either 
sample.  It  is  of  interest  to  note  that  NGC  4418  is  believed  to 
harbor  a  buried  AGN  based  on  its  mid-IR  spectral  proper¬ 
ties  (Roche  et  al.  1986;  Dudley  &  Wynn-Williams  1997; 
Spoon  et  al.  2001). 


We  note  that  our  convention  as  given  in  Table  2  for  calcu¬ 
lating  the  FIR  flux,  chosen  for  consistency  of  comparison 
with  previous  extragalactic  studies  of  the  [C  ii]  to  FIR  emis¬ 
sion  (e.g.,  Stacey  et  al.  1991;  Lord  et  al.  1996,  Malhotra 
et  al.  1997),  overestimates  the  actual  value  of  F[Cn|/+FIR, 
since  the  FIR  flux  estimate  accounts  for  the  FIR  emission 
only  between  42.5  and  122.5  /on  (Flelou  et  al.  1988).  The 
correction  factor  required  to  account  for  the  FIR  flux 
longward  of  122.5  /an  is  typically  in  the  range  of  1.4— 1.8 
(Sanders  &  Mirabel  1996).  Thus,  the  actual  [C  ii]  to  FIR 
ratio  is  slightly  less  than  presented  in  Figure  1.  For  some  of 
the  extended  normal  and  starburst  galaxies  for  which  the 
IRAS  Faint  Source  Catalog  (FSC)  may  have  missed  some 
of  the  FIR  emission  (Rice  et  al.  1988),  we  measured  the  ISO 
LWS  fluxes  using  ISAP,  again  from  42.5  to  122.5  //m  to 
match  the  IRAS  waveband,  and  plotted  these  values  in  Fig¬ 
ure  1.  For  some  of  the  more  extended  normal  and  starburst 
galaxies,  we  find  the  ISO  fluxes  to  be  as  much  as  40%  higher 
than  the  IRAS  FSC  fluxes.  For  galaxies  that  are  small  with 
respect  to  the  LWS  beam,  namely  all  ULIRGs  and  also  Arp 
299,  NGC  6240,  NGC  7552,  NGC  5713,  and  NGC  3256,  we 
used  the  IRAS  FSC  for  the  FIR  flux.  For  these  ULIRGs 
and  compact  non-ULIRGs,  the  ISO  LWS  and  IRAS  FSC 
flux  measurements  show  good  agreement  with  an  average 
ISO  to  IRAS  FIR  (42.5-122.5  /an)  flux  ratio  of  1.10  ±  0.09. 

For  comparison,  we  plot  in  Figure  2  the  6.2  /on  PAH 
feature  flux  (C.  C.  Dudley  et  al.,  in  preparation)  versus  the 
[C  ii]  emission.  The  galaxies  plotted  are  the  same  as  in  Fig¬ 
ure  1,  with  the  exception  of  IRAS  22491  —  1808,  which  is  not 
plotted.  Neither  [C  ii]  nor  the  6.2  /an  PAH  feature  were 
detected  for  this  ULIRG  and  therefore  its  upper  limits  do 
not  constrain  its  position  relative  to  the  band  shown  ( dashed 
lines).  The  6.2  /mi  emission  feature  arises  from  ionized 
PAHs  (Allamandola,  Hudgins,  &  Sandford  1999).  We  take 
the  6.2  /an  PAH  data  from  C.  C.  Dudley  et  al.  (in  prepara¬ 
tion),  who  measure  the  PAH  emission  from  ISO  PHOT-S 
spectra  and  ISO  CAM  CVF  spectra  integrated  over  the 
LWS  profile- weighted  beam  for  a  flux-limited  sample  of 
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Fig.  2. — 6.2  /an  PAH  feature  vs.  [C  ii]  158  /an  for  ULIRGs  (filled 
symbols)  and  normal/starburst  galaxies  (open  symbols).  The  galaxy  sample 
is  consistent  with  that  in  Fig.  1  (see  text  for  details).  The  6.2  /an  PAH  emis¬ 
sion  was  measured  from  ISO  CAM  CVF  spectra  ( open  symbols),  integrated 
over  the  full  LWS  beam  profile  for  galaxies  with  extended  emission,  and 
from  ISO  PHOT-S  spectra  (filled  symbols)  for  compact  galaxies  (C.  C. 
Dudley  et  al.,  in  preparation).  Statistical  error  bars  for  the  ISO  CAM  data 
are  shown,  although  in  virtually  all  cases  the  error  bars  are  smaller  than  the 
data  symbols.  Dashed  lines  show  a  range  of  a  factor  of  10. 
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Fig.  3. — 6.2  /an  PAH  to  [C  ii]  158  /im  line  ratio  distribution  for  the 
galaxies  from  Fig.  2.  Upper  limits  (3  a)  on  the  ratio  for  two  ULIRGs  are 
plotted  separately  and  marked  with  arrows.  The  upper  limit  for  IRAS 
20100—4156  is  26,  well  outside  the  plot  boundaries  and  not  a  useful  con¬ 
straint  on  the  ratio.  Also,  we  do  not  plot  galaxies  from  Fig.  1  with  AGN 
signatures  (NGC  6240,  NGC  4945,  NGC  1068,  and  Circinus),  in  order  to 
compare  the  ULIRG  distribution  with  that  of  normal  and  starburst 
galaxies  only. 

luminous  1R  galaxies.  For  three  of  the  extended  galaxies 
(Arp  299,  NGC  253,  and  Circinus),  only  ISO  CAM  data 
over  a  48"  x  48"  image  are  available  (smaller  than  the  69" 
LWS  aperture),  but  we  infer  from  larger  broadband  ISO 
CAM  images  that  these  images  include  most  of  the  flux. 
Dudley  et  al.  chose  the  6.2  pm  feature  because  it  is  relatively 
isolated  from  the  nearby  7. 7/8. 7  /mi  PAH  feature  blend, 
whose  measured  strength  can  be  affected  by  the  9.7  /an  sili¬ 
cate  absorption  feature.  Estimates  of  the  7.7  pi n  feature  flux 
must  therefore  account  for  the  effects  of  silicate  absorption, 
which  is  especially  difficult  for  redshifted  sources  with 
PHOT-S,  which  only  extends  to  11.5  /an  in  the  observed 
frame. 

In  Figure  3,  we  display  the  distribution  of  the  6.2  /an 
PAH  to  [C  n]  158  /an  line  ratio  in  histogram  format  for  the 
ULIRGs  and  for  the  normal  and  starburst  galaxies  (no  gal¬ 
axies  with  AGN  signatures  are  included  in  the  comparison 
sample).  Both  Figures  2  and  3  show  that,  in  contrast  to  the 
[C  ii]  to  FIR  deficit  found  in  ULIRGs,  the  PAH  to  [C  ii] 
ratio  is  similar  in  ULIRGs  and  the  comparison  sample,  i.e., 
both  [C  ii]  and  the  6.2  /an  PAH  feature  display  a  propor¬ 
tionally  similar  deficit  with  respect  to  FIR  in  ULIRGs  and 
normal/ starburst  galaxies.  Using  PHOT-S,  Lutz  et  al. 
(1998)  and  Rigopoulou  et  al.  (1999)  find  the  7.7  /an  PAH  to 
IRAS  60  /mi  ratio  emission  to  be  approximately  3  times 
weaker  in  ULIRGs  than  in  starbursts.  For  these  studies,  the 
starburst  galaxies  are  mostly  relatively  nearby  sources, 
where  in  Rigopoulou  et  al.,  11  out  of  14  starburst  galaxies 
are  within  cz  560  km  s  1 .  Since  the  aperture  of  the  PHOT- 
S  observations  was  24"  x  24",  while  the  IRAS  60  /mi  fluxes 
likely  arise  in  a  significantly  larger  area  in  the  starburst  gal¬ 
axies,  the  factor  of  3  difference  in  the  PAH  to  IRAS  60  pm 
ratio  is  likely  an  underestimate.  Here  we  attempted  to  avoid 
this  problem  by  using  ISO  CAM  CVF  measurements  for 
the  nearby  starburst  galaxies,  from  which  it  is  possible  to 
measure  only  the  PAH  emission  corresponding  to  the  LWS 
measurement.  Furthermore,  for  the  reasons  stated  above, 
we  used  the  6.2  pm  feature  rather  than  the  7.7  /mi  feature 
used  by  Rigopoulou  et  al.  Thus,  we  believe  that  the  similar¬ 
ity  between  ULIRGs  and  starbursts  in  the  PAH  to  [C  ii] 
ratios  shown  here  is  unbiased  in  these  important  respects. 


3.  DISCUSSION 

In  Paper  I,  we  discussed  three  possible  explanations  for 
the  [C  ii]  deficit:  (1)  self-absorbed  or  optically  thick  [C  n] 
emission;  (2)  saturation  of  the  [C  n]  emission  in  a  photodis- 
sociated  gas  with  a  high  ratio  of  incident  UV  flux  G0  to  n;  or 
(3)  the  presence  of  a  soft  UV  radiation  field  suggesting,  for 
example,  an  older  stellar  population  or  absorption  of  UV  in 
dusty  photoionization  regions.  Based  on  the  [C  n]  observa¬ 
tions  alone,  we  could  not  make  a  definitive  case  in  favor  of 
only  one  of  the  three  proposed  explanations,  although  self- 
absorbed  or  optically  thick  [C  n]  emission  seemed  the  least 
likely  explanation  for  the  [C  ii]  deficit.  Using  the  added  ISO 
[C  ii]  and  [O  i]  observations,  we  explore  these  and  other 
possible  sources  for  the  [C  n]  deficit  not  posed  in  Paper  I. 

3.1.  Saturation  of[ C  ii]/  High  6/ / n  PDRs 

For  PDRs  with  a  high  ratio  of  Go  to  n  (Go/«M-10  cm3), 
theoretical  PDR  models  show  that  [C  ii]  saturates  relative  to 
FIR  such  that  the  ratio  of  Fjc„]  to  Ffir  drops  below  values  of 
1CU4  (e.g.,  Wolfire,  Tielens,  &  Hollenbach  1990;  Kaufman 
et  al.  1999).  Here,  the  high  G0  increases  the  positive  charge  of 
dust  grains,  thereby  reducing  the  energy  of  ejected  photo¬ 
electrons  and  the  fraction  of  UV  photon  energy  that  is  con¬ 
verted  to  gas  heating.  As  a  result,  gas  cooling  via  the  [C  n]  line 
falls.  (We  note  that  another  type  of  [C  ii]  saturation  occurs 
when  the  gas  density  n  lies  well  above  the  [C  ii]  critical 
density.  We  discuss  this  case  further  in  §  3.2.1.) 

Our  observed  [C  ii]  line  fluxes  combined  with  the 
available  [O  i]  observations  allow  us  to  measure  n  and  Go  for 
several  sample  ULIRGs  and  thus  directly  address  the  possi¬ 
bility  of  saturation  of  the  [C  n]  emission  in  PDR  gas.  The 
contour  lines  shown  in  Figure  4  depict  the  expected 
behavior  of  [C  ii]  158  //m/[0  i]  63  pm  and  ([C  ii]  1 58  /mi  + 
[O  i]  63  /im)//FIR  as  a  function  of  n  and  G0  for  extragalactic 
PDRs.  (The  UV  flux  Go  is  normalized  to  the  solar  neighbor¬ 
hood  radiation  field,  or  1.6  x  10-3  ergs  cm~2  s_1;  Habing 
1968.)  The  contour  plot  is  based  on  the  predictions  of  the 
Kaufman  et  al.  (1999)  PDR  code,  which  is  an  updated  ver¬ 
sion  of  earlier  works  stemming  from  the  study  by  Tielens  & 
Hollenbach  (1985).  Among  other  updates  to  the  input 
physics,  the  code  includes  the  effects  of  gas  heating  associ¬ 
ated  with  the  photoelectric  ejection  of  electrons  from  very 
small  grains  (<15  A)  and  PAHs  (Bakes  &  Tielens  1994). 
The  standard  PDR  model  of  Kaufman  et  al.  adopts  the  frac¬ 
tional  gas  phase  abundances  as  measured  by  Sembach  & 
Savage  (1996;  i.e.,  C/H  =  1.4  x  1(U4,  O/H  =  3  x  1()-4). 
For  an  extragalactic  PDR,  approximately  half  of  the 
observed  IR  dust  continuum  emission  comes  from  clouds 
facing  away  from  the  observer  that  do  not  contribute  to  the 
[C  ii]  158  pm  and  [O  i]  63  pm  emission  directed  along  the  line 
of  sight  owing  to  optical  depth  effects  (see  Wolfire  et  al.  1990 
and  the  Appendix  of  this  paper).  This  “  excess  ”  IR  contin¬ 
uum  emission  is  accounted  for  in  Figure  4  (and  Fig.  5)  by 
setting  /FIR  equal  to  twice  the  theoretical  continuum  inten¬ 
sity  emitted  by  a  face-on  PDR  [or  /FjR  =  2  x  (2  x  1.6  x 
10~3G0/47t)  =  5.1  x  10_4G0  ergs  cnu2  s-1  sr-1]. 

In  order  to  use  Figure  4  to  measure  n  and  Go  for  our 
sample  of  ULIRGs,  the  observed  [C  n],  [O  i],  and  FIR 
emission  must  arise  predominantly  from  the  observed 
ensemble  of  PDR  surfaces.  (Since  the  physical  size  of  the 
beam  at  the  distance  to  the  source  is  large  compared  to  a 
single  PDR,  the  observed  [C  ii],  [O  i],  and  FIR  emission 
arises  from  many  PDR  surfaces  separated  in  both 
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Fig.  4. — Plot  of  the  [C  ii]  158  //m/[0 1]  63  /im  ratio  vs.  the  ratio  of  ([C  n]  158  /im  +  [O  i]  63  /im)//r[R  emitted  from  the  surface  of  a  PDR.  Models  of  Kaufman 
et  al.  (1999)  are  used  with  the  standard  model  parameter  set  and  assuming  that  /fir  =  2  x  (2  x  1.3  x  1(D4)Go  ergs  cm-2  s-1  sr-1.  The  five  ULIRGs — 
Mrk  231,  Mrk  273,  IRAS  17208-0014,  IRAS  20551-4250,  and  IRAS  23365+3604— from  this  study  with  both  [C  n]  158  //m  and  [O  i]  63  /un  emission  line 
measurements  are  shown  (filled  symbols).  Derived  propagation  errors  resulting  front  errors  in  the  line  and  FIR  flux  measurements  imply  variations  in  the 
depicted  values  of  n  and  Go  of  no  more  than  a  factor  of  a  few.  For  comparison,  we  plot  the  objects  (open  circles)  in  the  Malhotra  et  al.  (2001)  sample  of  60 
normal  and  starburst  galaxies.  For  one  of  the  points  of  the  Malhotra  et  al.  sample  we  depict  sample  error  bars,  which  are  constant  for  all  points.  As  with 
the  ULIRG  sample,  we  plot  only  those  objects  front  the  Malhotra  et  al.  sample  for  which  both  [C  n]  158  /im  and  [O  i]  63  /mi  lines  were  detected.  For  Mrk  231, 
the  vector  depicts  the  shift  in  the  location  of  this  galaxy  in  the  figure  if  we  correct  for  the  maximum  allowable  decrease  in  the  [O  i]  63  jum  emission  owing  to 
self-absorption. 


physical  and  velocity  space.)  Some  portion  of  the  diag¬ 
nostic  line  emission  may  arise  from  non-PDR  gas.  (For 
now,  we  neglect  non-PDR  sources  for  the  FIR  emission; 
we  treat  this  case  further  in  §  3.2.)  For  example,  ionized 
gas  in  the  form  of  classical  H  ii  regions  or  warm 
(T  ~  8  x  103  K)  diffuse  gas  can  emit  [C  ii]  radiation;  like¬ 
wise,  shock-excited  gas  can  generate  significant  [O  i] 
emission.  The  fraction  of  non-PDR  line  emission  should 
be  subtracted  before  using  Figure  4. 

Although  we  are  unable  to  estimate  the  fractional  contri¬ 
bution  from  non-PDR  gas  for  the  sources  studied,  we  argue 
that  the  contribution  is  less  than  that  of  the  PDR  compo¬ 
nent.  One  method  for  estimating  the  [C  ii]  contribution  from 
ionized  gas  is  to  rely  on  H  n  region  modeling,  as  presented 
by  Carral  et  al.  (1994)  and  Satyapal  et  al.  (2003).  Carral 
et  al.  point  out  that  extragalactic  H  ii  regions  can  contribute 
as  much  as  30%-50%  of  the  total  [C  n]  luminosity.  Using  the 
[N  ii]  122  /on  line,  which  traces  ionized  gas  in  both  H  ii 
regions  and  diffuse  media,  Malhotra  et  al.  (2001)  and 
Negishi  et  al.  (2001)  also  attribute  as  much  as  30%-50%  of 
the  [C  ii]  emission  to  ionized  gas.  For  the  UFIRGs  studied 
here,  [N  n]  122  //m  measurements  are  not  available.  In  the 
case  of  the  prototypical  luminous  starburst  galaxy  Arp  299 
(Lir  ~  4.8  x  1011  L0),  with  luminosity  approaching  that  of 
the  UFIRGs,  Satyapal  et  al.  find  that  the  contribution  from 
ionized  gas  to  the  [C  ii]  luminosity  is  closer  to  1 5%  based  on 
the  relative  [C  ii]  and  [N  ii]  122  /on  line  strengths  in  Arp  299. 
Satyapal  et  al.  conclude  that,  compared  to  the  warm  diffuse 
gas,  classical  H  n  regions  provide  the  bulk  of  the  ionized  gas 
contribution  to  the  [C  ii]  emission. 

For  the  shock  contribution  to  the  [O  i]  emission,  shock 
models  (e.g.,  Hollenbach  &  McKee  1989)  predict  [O  i] 


63  /on  to  [C  ii]  158  /on  ratios  greater  than  3.  Such  ratios  are 
roughly  consistent  with  those  seen  in  AGNs  and  starbursts 
with  a  strong  shock-excited  component  (e.g.,  Spinoglio  & 
Malkan  1992;  van  der  Werf  et  al.  1993).  Based  on  the  gener¬ 
ally  low  ratio  (~1)  of  [O  i]  to  [C  n]  emission  in  the  UFIRGs, 
however,  we  believe  that  the  contribution  to  the  [O  i] 
emission  from  shock-excited  gas  is  small. 

In  Figure  4,  we  plot  five  UFIRGs — Mrk  231,  Mrk  273, 
IRAS  17208—0014,  IRAS  20551-4250,  and  IRAS 
23356+3604 — for  which  we  have  both  [C  ii]  158  //m  and 
[O  i]  63  /mi  emission  line  measurements.  For  comparison, 
we  also  plot  the  objects  in  the  Malhotra  et  al.  (2001)  sample 
of  60  normal  and  starburst  galaxies.  With  the  possible 
exception  of  Mrk  23 1 ,  we  see  that  the  implied  PDR  gas  den¬ 
sity  n  (~100-103  emu3)  is  similar  to  that  of  the  normal  and 
starburst  galaxies,  i.e.,  the  [C  ii]  to  [O  i]  63  /mi  ratios  for  the 
ULIRGs  are  comparable  to  those  of  the  normal  and  star- 
burst  galaxies  shown.  Although  the  derived  PDR  densities 
are  comparable  to  the  normal  and  starburst  galaxies,  they 
fall  well  below  the  interstellar  gas  densities  implied  by  CO 
observations  of  the  ULIRGs  (Solomon  et  al.  1997). 
Whereas  the  typical  PDR  density  of  the  ULIRGs  plotted  in 
Figure  4  is  ~200  cm~3,  the  interstellar  gas  densities  mea¬ 
sured  by  Solomon  et  al.  range  from  ~700  emu3  (Mrk  231 
and  IRAS  23365+3604)  to  1.3-1. 6  x  103  cm-3  (Mrk  273 
and  IRAS  17208-0014). 

In  Figure  4,  the  implied  UV  flux  Go  (~500-103)  incident 
on  the  PDR  gas  is  significantly  higher  than  that  of  the  nor¬ 
mal  and  starburst  galaxies  by  a  factor  of  2-5.  As  mentioned, 
for  PDRs  with  high  Go  to  n  ratios,  the  high  UV  flux 
increases  grain  charging  and  reduces  the  efficiency  of  photo¬ 
electric  heating,  and  thus  [C  ii]/FIR  falls  (Wolfire,  Tielens, 
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&  Hollenbach  1990;  Kaufman  et  al.  1999).  Correcting  the 
[C  ii]  and  [O  i]  line  fluxes  for  contributions  from  non-PDR 
sources  to  the  [C  ii]  and  [O  i]  emission  could  alter  the  derived 
n  and  G0  considerably,  but  we  would  not  expect  the  G0  to  n 
ratio  to  decrease.  For  example,  applying  a  reasonable  cor¬ 
rection  to  the  [C  ii]  flux  to  account  for  emission  from  ionized 
gas  would  move  the  ULIRG  data  points  in  Figure  4  down 
and  to  the  left  in  such  a  way  that  both  n  and  Go  would 
increase  by  a  similar  amount.  Fikewise,  lowering  the  [O  i] 
flux  by  some  correction  factor  to  account  for  the  shocked 
gas  would  move  the  UFIRG  data  points  toward  higher  G0 
and  lower  n,  thus  potentially  increasing  Go,//;.  For  large  cor¬ 
rections  in  the  [O  i]  emission  of  a  factor  of  2  or  more,  the 
derived  G(,/«.  which  ranges  from  1  to  5  in  Figure  4,  would 
increase  by  a  factor  of  a  few. 

A  high  G0/  n  cause  for  the  UFIRG  [C  ii]  deficit  would 
agree  with  the  conclusion  of  Malhotra  et  al.  (1997,  2001), 
who  suggest  that  high  Go/  n  PDRs  (specifically,  high  n  and 
higher  Go)  account  for  the  [C  ii]  flux  deficiencies  seen  in  the 
FIR-bright  objects  among  their  sample  of  normal  galaxies. 
Malhotra  et  al.  base  their  argument  on  the  observed 
decrease  in  the  [C  ii]  to  FIR  ratio  as  a  function  of  the  60-100 
//in  flux.  The  60-100  //nr  emission  from  PDRs  measures 
average  temperatures  of  the  large  grain  population  and 
should  therefore  increase  with  increasing  UV  flux  Go-  In 
Figure  5,  we  plot  the  Figure  1  data  as  a  function  of  the  60- 
100  /zm  flux  ratio.  Here,  the  drop  in  [C  ii]  to  FIR  flux  with 
the  60-100  //m  flux  ratio  shows  that  the  full  UFIRG  sample 
data  are  consistent  with  the  result  found  by  Malhotra  et  al. 
We  emphasize,  however,  that  it  is  a  decrease  in  [C  n]/FIR 
with  increasing  G0/«,  not  with  increasing  G0  alone,  that 
points  to  a  decrease  in  the  photoelectric  heating  efficiency  as 
the  cause  for  the  [C  ii]  deficit.  As  pointed  out  by  Negishi 
et  al.  (2001),  we  would  not  expect  [C  n]/FIR  to  decrease  if 
both  G0  and  n  increase  together.  If  we  take  Figure  4  at  face 
value,  then  G0  would  seem  to  increase  with  respect  to  n  in 
UFIRGs  as  compared  to  normal  and  starburst  galaxies; 
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Fig.  5. — [C  ii]  158  // m  line  to  FIR  flux  ratio  vs.  the  IRAS  60-100  //m  flux 
ratio  for  all  galaxies  plotted  in  Fig.  1 .  The  filled  circles  denote  the  ULIRGs, 
the  squares  denote  the  comparison  normal  and  starburst  galaxies,  and  the 
gray-shaded  squares  denote  the  galaxies  among  the  comparison  sample 
with  AGN  signatures.  The  error  bars  for  the  ULIRGs  are  statistical;  for  the 
comparison  galaxies,  the  statistical  errors  bars  are  smaller  than  the  symbols 
and  are  therefore  not  plotted.  The  arrows  represent  3  a  upper  limits.  For 
extended  galaxies,  the  60-100  //m  ratios  were  derived  from  the  Rice  et  al. 
(1988)  IRAS  catalog  rather  than  the  IRAS  Faint  Source  Catalog.  For 
galaxies  with  no  IRAS  data  (e.g.,  Circinus),  the  60-100  //m  flux  ratio  was 
derived  from  the  LWS  full  spectrum  using  the  photometry  routine  in  ISAP. 


thus,  Figure  5  appears  to  be  consistent  with  a  high  G()/n 
origin  for  the  [C  ii]  deficiency. 

Assuming  that  we  are  correct  in  applying  the  PDR 
models  to  the  UFIRGs,  we  explore  the  degree  to  which 
uncorrected  self-absorption  in  the  [O  i]  63  /;m  line  could 
influence  the  location  of  the  UFIRGs  in  Figure  4.  Fischer 
et  al.  (1999)  have  reported  such  self-absorption  in  the  [O  i] 
63  /zm  line  in  the  prototypical  UFIRG  Arp  220.  The  upper 
limit  on  the  [O  i]  145  /zm  line  places  a  strong  constraint  on 
the  amount  by  which  we  can  feasibly  correct  the  [O  i]  63  //nr 
line  for  self-absorption.  Unlike  the  ground-state  [O  i]  63  //nr 
line,  which  is  seen  in  absorption  in  Arp  220,  the  [O  i]  145  //nr 
line  should  not  be  self-absorbed,  since  it  originates  from  an 
energy  level  that  is  228  K  above  the  ground  state  and  should 
not  be  populated  in  cold  or  diffuse  gas.  If,  for  example,  we 
increase  the  [O  i]  63  /zm  line  emission  by  a  factor  of  2, 
Mrk  231  shifts  to  a  location  in  Figure  4  as  shown  by  the 
absorption  vector  corresponding  to  n  ~  700  cnu3  and 
Go  ~  3  x  103.  According  to  Figure  6,  for  this  value  of  n  and 
G0,  the  [C  ii]  158  //m  to  [O  i]  145  /zm  ratio  is  ~7,  which  corre¬ 
sponds  to  the  measured  lower  limit  on  the  [C  ii]  to  [O  i]  145 
/zm  ratio  for  Mrk  23 1 .  Any  correction  for  [O  i]  63  /zm  self¬ 
absorption  greater  than  a  factor  of  2  would  therefore  yield 
a  predicted  [C  ii]  to  [O  i]  145  //nr  smaller  than  ~7,  in 
disagreement  with  the  observed  lower  limit  for  this  ratio. 

When  we  apply  this  same  constraint  to  the  other  UFIRGs 
in  Figure  4  with  [O  i]  145  /zm  upper  limits,  we  find  that  even 
small  corrections  for  self-absorption  lead  to  inconsistencies 
with  the  observations,  i.e.,  the  [O  i]  63  /zm  line  is  likely  not 
self-absorbed  in  these  UFIRGs.  For  Mrk  273,  the  observed 
lower  limit  on  the  [C  n]  to  [O  i]  145  //m  ratio  is  18.  According 
to  Figure  6,  a  [C  ii]  to  [O  i]  145  /zm  ratio  of  18  matches  the 
derived  n  (~200  cnr3)  and  6/  (~700)  from  Figure  4.  Any 


Fig.  6. — Ratio  of  the  [C  n]  158  n m  line  intensity  to  the  [O  i]  145  //m  line 
intensity  emitted  from  the  surface  of  a  PDR  as  a  function  of  the  cloud 
density  n  and  the  UV  flux  incident  on  the  cloud  Go.  Models  of  Kaufman 
et  al.  ( 1 999)  are  used  with  the  standard  model  parameter  set. 
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upward  correction  to  the  [O  i]  63  /an  line  flux  would  drive 
up  n  and  G0,  thus  dropping  the  predicted  [C  ii]  to  [O  i] 
145  /an  line  ratio  below  the  measured  lower  limit,  again  in 
disagreement  with  the  observations.  The  same  holds  true 
for  IRAS  17208—0014,  the  remaining  ULIRG  in  Figure  4 
with  an  observed  [O  i]  145  /on  upper  limit.  Thus,  with  the 
possible  exception  of  Mrk  231,  [O  i]  self-absorption  does 
not  appear  to  affect  the  derived  PDR  gas  densities  and  UV 
fluxes  for  the  ULIRGs  in  Figure  4. 

Irrespective  of  the  potential  effects  of  [O  i]  63  /on  self¬ 
absorption  on  the  derived  PDR  parameters,  the  [O  i] 
145  /mi  line,  together  with  the  [C  ii]  and  FIR  fluxes,  con¬ 
strains  the  derived  n  and  60.  For  example,  for  Arp  220, 
T[cii]/T[oi]i45,m>ll,  T"[c „] /Tfir  =  2.09  x  lCU4;  for  Mrk 
273,  F\cu]/Fyo ,]  i45/mi>18,  Cpuj/CpiR  =  5.37  x  KU4;  and  for 
IRAS  17208-0014,  F[Cu]/F[0l]l45fim^U,  F[Cu]/Ffir  = 
4.73  x  10  4.  The  only  PDR  regimes  with  these  line  ratios 
that  are  in  keeping  with  both  Figures  6  and  7  are  either  (1) 
Go  -  (0.7-2)  x  103  and  n  -  100  cm-3  or  (2)  G0  -  1-10  and 
n  —  105  cm-3.  The  former  case  would  characterize  a  high 
Go/h  PDR,  but  n  —  100  cm-3  is  much  lower  than  the 
derived  average  interstellar  densities  of  — 104  cm-3  for  Arp 
220,  -1300  cm-3  for  Mrk  273,  and  -1600  cm-3  for  IRAS 
17208—0014  from  CO  observations  (Solomon  et  al.  1997; 
Scoville,  Yun,  &  Bryant  1997).  (As  we  have  already  noted, 
the  observations  and  PDR  models  together  significantly 
underestimate  the  PDR  gas  density  for  the  ULIRGs  in  Fig. 
4  as  compared  to  previous  CO  estimates  of  the  interstellar 
gas  density.)  The  density  for  the  latter  case  is  considerably 
higher,  on  the  other  hand,  but  values  of  G0  —  1-10  and 
G0/n  -  10-4  in  Arp  220,  Mrk  273,  and  IRAS  17208-0014 
are  highly  unlikely  given  the  high  luminosity  of  these  gal- 


Fig.  7. — Ratio  of  the  [C  ii]  158  fim  line  intensity  to  the  FIR  continuum 
intensity  emitted  from  the  ensemble  of  extragalactic  PDRs  as  a  function  of 
the  cloud  density  n  and  the  UV  flux  incident  on  the  cloud  Go.  Models  of 
Kaufman  et  al.  (1999)  are  used  with  the  standard  model  parameter  set. 


axies.  Thus,  we  find  that  any  reasonable  PDR  model  that 
can  account  for  the  observed  [C  n],  [O  i]  63  and  145  /mi, 
and  FIR  emission  does  not  compare  well  with  previous 
observations  of  the  neutral/molecular  ISM  in  ULIRGs. 

Putting  aside  for  the  moment  the  disparity  between  the  n 
and  Go  for  ULIRGs  derived  from  the  PDR  analysis  above 
and  from  previous  studies,  we  examine  the  issue  of  high 
Go/n  PDRs  and  their  effect,  if  any,  on  the  observed  similar¬ 
ity  of  the  6.2  /on  PAH  to  [C  ii]  ratio  in  ULIRGs  and  in  nor¬ 
mal  and  starburst  galaxies  seen  in  Figures  2  and  3.  In  order 
for  high  Go/n  PDRs  to  characterize  ULIRGs  but  not  most 
normal/ starburst  galaxies,  high  Go/n  PDRs  must  decrease 
both  the  [C  ii]  and  PAH  flux  together  relative  to  FIR  in 
ULIRGs,  excluding  for  now  other  effects  such  as  mid-IR 
extinction  (discussed  in  §  3.4),  which  could  act  to  reduce  the 
PAH  flux  relative  to  [C  ii].  In  Figure  8,  we  plot  the  ratio  of 
the  6.2  /mi  PAH  to  [C  ii]  flux  versus  the  60-100  /mi  ratio, 
thought  to  scale  with  Go  in  the  PDR.  (As  with  Fig.  2,  IRAS 
22491  —  1808  is  not  included  in  this  plot  since  neither  its  6.2 
/mi  PAH  nor  [C  ii]  line  emission  have  been  detected.)  We  see 
that  the  PAH  to  [C  n]  ratio  stays  fairly  constant  over  the 
entire  range  of  60-100  /mi  ratios  characteristic  of  normal/ 
starburst  galaxies  and  UFIRGs.  We  note  that  two  of  the 
four  AGN-signature  galaxies  are  relatively  low  in  6.2  /mi 
PAH  to  [C  ii]  flux;  this  tendency  in  AGNs,  if  real,  may  be 
due  to  an  enhanced  contribution  from  a  non-PAH  generat¬ 
ing  component  (such  as  the  narrow-line  region),  AGN  proc¬ 
essing  of  PAH  molecules,  or  other  AGN  effects  (Imanishi  & 
Dudley  2000).  Helou  et  al.  (2001)  also  observed  an  implied 
constancy  in  the  PAH  to  [C  ii]  ratio  for  their  ISO  normal 
galaxies  key  project  using  broadband  mid-IR  photometry 
to  measure  the  PAH  emission.  (We  note  that  low-metallicity 
dwarfs  can  display  systematically  weaker  PAH  to  [C  ii] 
ratios;  Madden  2000;  Hunter  et  al.  2001.)  Helou  et  al. 
(2001)  explain  the  close  association  of  PAH/[C  ii]  with  60- 
100  /mi  color  by  assuming  that  the  PAH  emission  measures 
the  PDR  photoelectric  heating,  and  [C  ii]  measures  the  gas 
cooling.  In  this  scenario,  the  decrease  in  [C  n]/FIR  (Fig.  5) 
and  PAH/FIR  (inferred  from  Figs.  5  and  8)  with  increasing 
radiation  field  would  suggest  a  decrease  in  the  importance 
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Fig.  8. — 6.2  ;/m  PAH  to  [C  ii]  158  fim  line  flux  ratio  vs.  the  IRAS  60-100 
/im  flux  ratio  for  all  galaxies  plotted  in  Fig.  2.  The  filled  circles  denote  the 
ULIRGs,  the  squares  denote  the  comparison  normal  and  starburst 
galaxies,  and  the  gray-shaded  squares  denote  the  galaxies  among  the  com¬ 
parison  sample  with  AGN  signatures.  The  error  bars  and  upper  limits  are  1 
and  3  cr  values,  respectively.  The  60-100  fim  ratios  are  derived  as  in  Fig.  5. 
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of  PAHs  and  small  grains  relative  to  large  grains  in  both  gas 
heating  and  dust  cooling. 

Although  [C  ii]  as  a  measure  of  PDR  gas  cooling  is 
well  established  via  both  theory  and  observations,  the 
association  between  the  6.2  pm  feature  emission  from 
ionized  PAHs  and  the  photoelectric  process  as  a  function 
of  G0  is  not  so  definitive.  High  G0/«  conditions  causing 
low  values  of  the  [C  ii]  to  FIR  ratio  would  not  produce  a 
correspondingly  low  6.2  pm  feature  flux,  unless  a  signifi¬ 
cant  fraction  of  the  grains  responsible  for  the  6.2  /on 
feature  are  destroyed  by  these  conditions.  However,  labo¬ 
ratory  studies  indicate  that  PAH  sizes  larger  than  50  car¬ 
bon  atoms  probably  dominate  the  emission  in  this  band 
(Hudgins  &  Allamandola  1999),  while  theoretical  studies 
indicate  that  it  is  difficult  to  destroy  PAHs  with  more 
than  50  carbon  atoms  in  the  regimes  of  G^/n  derived  here 
(Allain,  Leach,  &  Sedlmayr  1996a,  1996b).  According  to 
Bakes  &  Tielens  (1994),  it  would  appear  that  grains  up  to 
at  least  100  A  in  size  would  have  to  be  destroyed  in  order 
to  lower  the  heating  by  a  factor  of  10.  Recent  ISO  obser¬ 
vations  also  show  little  change  in  the  mid-IR  PAH  spec¬ 
tral  shapes  over  a  wide  range  of  G0  ( 1  10s;  Uchida  et  al. 
2000),  even  though  laboratory  studies  show  significant 
changes  in  relative  strengths  of  different  features  from 
ionized  to  neutral  PAHs.  Because  of  such  results,  it  is  dif¬ 
ficult  to  unambiguously  associate  the  6.2  pm  feature 
emission  from  ionized  PAHs  with  the  photoelectric  proc¬ 
ess  as  a  function  of  G0,  even  though  Figure  8  and  the 
similar  finding  of  Helou  et  al.  suggest  that  PAHs  closely 
track  the  photoelectric  heating. 

Before  we  conclude  our  discussion  of  high  Go/n  PDRs 
and  the  relative  PAH  and  [C  ii]  emission,  we  note  another 
interesting  result  of  the  Helou  et  al.  (2001)  study.  Helou 
et  al.  find  that  at  larger  values  of  the  60-100  /mi  flux  ratio, 
the  [O  i]  63  pm  emission  overtakes  [C  ii]  in  cooling  lumin¬ 
osity.  Helou  et  al.  point  out  that  this  observation  conflicts 
with  the  argument  that  PAH/[C  ii]  is  constant  because  it  is 
tied  to  the  photoelectric  effect,  since  as  G0  increases,  the 
ratio  of  ([C  ii]  +  [O  i])/PAH  rises,  suggesting  an  increase  in 
the  photoelectric  heating  efficiency  of  the  PAHs,  which  con¬ 
tradicts  the  physics.  One  scenario  devised  by  Helou  et  al.  to 
explain  the  observations  is  a  two-component  model 
whereby  the  FIR  emission  conies  from  both  an  “active” 
star  formation  PDR  component  and  a  “quiescent”  PDR 
component.  The  “  quiescent  ”  PDR  component  would  dom¬ 
inate  at  low  values  of  60  /on/ 100  //m  and  display  high  values 
of  [C  n]/[0  i],  PAH/FIR,  and  [C  ii]/FIR;  the  “active” 
component  would  dominate  at  high  60  //m/100  p m  and 
show  low  values  of  [C  n]/[0  i],  PAH/FIR,  and  [C  ii]/F1R. 
Although  we  do  not  observe  especially  low  values  of  [C  n]/ 
[O  i],  a  two-component  model  for  the  FIR  emission  raises 
the  possibility  that  the  ratio  of  60  //m/100  /mi  as  plotted  in 
Figures  5  and  8  may  not  necessarily  track  Go  in  the  region 
that  gives  rise  to  the  bulk  of  the  PAH  and  [C  ii]  emission. 
Instead,  another  FIR  component  could  contribute  to  the 
integrated  60  and  100  /zm  fluxes,  producing  warmer  60  pm/ 
100  fim  colors  and  increasing  the  FIR  flux  relative  to  [C  ii] 
and  PAH  as  observed  in  Figures  5  and  8.  The  existence  of 
multiple  components  to  the  FIR  flux  would  bring  into  ques¬ 
tion  our  application  of  the  one-dimensional,  plane-parallel 
PDR  models  presented  in  this  section.  In  the  following  sec¬ 
tion,  we  explore  in  more  detail  the  possible  contribution  of 
non-PDR  components  to  the  FIR  emission  and  its  bearing 
on  the  [C  ii]  deficit. 


3.2.  Applicability  of  PDR  Models: 

Non-PDR  FIR  Components? 

A  key  assumption  underlies  the  one-dimensional,  plane- 
parallel  PDR  analysis  presented  in  §  3.1:  the  beam  area  fill¬ 
ing  factors  for  the  [C  n],  [O  i],  and  FIR  emission  are  taken  to 
be  equal  or  nearly  so.  Simply  put,  the  PDR  analysis  assumes 
that  the  bulk  of  the  [C  ii],  [O  i],  and  FIR  emission  comes 
from  the  same  region  of  PDR  gas  with  roughly  similar  aver¬ 
age  PDR  parameters  and  that  no  other  interstellar  compo¬ 
nents  contribute  significantly  to  the  observed  emission. 
Here,  the  beam  area  filling  factor  $  refers  to  the  ratio  of  the 
effective  PDR  surface  area  intercepted  by  the  beam  to  the 
beam  projected  area.  Multiple  PDRs  along  a  line  of  sight 
can  contribute  to  the  filling  factor  as  long  as  the  background 
emission  can  escape  to  enter  the  beam.  This  can  be  due  to 
velocity  shifts  between  PDR  layers  so  that  line  emission  is 
not  blocked  by  foreground  gas,  or  to  optically  thin  line 
emission  at  the  same  velocity.  For  distant  extragalactic 
sources  where  the  beam  encompasses  the  entire  galaxy,  fill¬ 
ing  factor  effects  may  be  important  (e.g.,  Wolfire, 
Hollenbach,  &  Tielens  1989).  If  non-PDR  components  of 
line  or  continuum  emission  are  present  in  ULIRGs  or  if 
multiple  or  varying  PDR  components  with  vastly  different 
conditions  are  present,  then  the  PDR  models  implying  high 
G0/n  PDRs  do  not  apply.  Such  scenarios  could  account  for 
the  measured  deficit  in  [C  n]  relative  to  FIR,  since  once  cor¬ 
rected  for  filling  factor  effects,  the  [C  ii]  to  FIR  ratios  in  the 
ULIRGs  may  resemble  the  ratios  in  normal  and  starburst 
galaxies. 

We  can  estimate  the  [C  ii]  beam  area  filling  factor 
<b[c  h]  in  the  four  ULIRGs  for  which  we  estimated  n  and 
Go  above  by  comparing  the  observed  [C  n]  intensity  to 
the  predicted  [C  ii]  intensity  for  the  derived  n  and  G0. 
Table  3  shows  the  observed  [C  ii]  intensities  for  Mrk  273, 
IRAS  17208-0014,  IRAS  20551-4250,  and  IRAS 
23365+3604.  We  derive  the  [C  ii]  intensities  from  the 
observed  fluxes  as  measured  in  a  69"  beam.  Along  with 
the  observed  intensities,  Table  3  gives  the  predicted  [C  ii] 
intensities  (Kaufman  et  al.  1999)  for  a  unity  beam  area 
filling  factor  given  n  and  G0  as  derived  in  the  previous 
section.  The  ratio  of  the  measured  [C  ii]  intensity  to  the 
predicted  [C  n]  intensity  then  provides  an  estimate  of  the 
beam  area  filling  factor  for  the  observed  [C  ii]  (Table  3, 
col.  [4]).  We  note  that  although  the  beam  area  filling  fac¬ 
tors  derived  in  this  manner  are  rough  estimates,  this 
approach  can  be  used  to  make  a  relative  comparison  of 
beam  area  filling  factors  for  various  emission  line  or 
continuum  components,  as  illustrated  in  the  following 
sections. 

3.2.1.  CO  ( 1-0)  Emission  Component 

Having  derived  4>[C  „]  for  four  ULIRGs  in  our  sample,  we 
now  examine  whether  the  filling  factors  of  other  PDR  emis¬ 
sion  lines,  and  ultimately  the  FIR  continuum,  match  those 
of  [C  ii].  This  filling  factor  analysis  is  in  effect  a  consistency 
check  on  our  assumptions  that  the  bulk  of  the  PDR  line 
diagnostics  and  FIR  emission  come  from  the  same  regions 
of  PDR  gas  with  roughly  similar  average  PDR  parameters. 
To  our  knowledge,  such  an  analysis  is  done  here  in  detail  for 
the  first  time.  We  begin  with  the  2.6  mm  CO  (1-0)  emission. 
Kaufman  et  al.  (1999)  model  the  behavior  of  the  CO  (1-0) 
intensity  as  a  function  of  n  and  G0,  again  assuming  unity 
beam  area  filling  factor.  Using  the  Kaufman  et  al. 
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TABLE  3 

[C  n]  Beam  Area  Filling  Factors  4>[C  nj  for  Four  ULIRGs  and  Comparison 
Normal/Starburst  Galaxies 


Observed  /[Cn]  158  Predicted  /[Cn]  i5g  for  %  „]  =  1 
Source  (ergs  cm-2  s-1sr-1)  (ergscm-2s_,sr-1)  (E>[cri] 


ULIRGs 

Mrk  273 .  6.0  x  10“6  2  x  10~4  0.030 

IRAS  17208-0014 .  7.8  x  10“6  2  x  10~4  0.039 

IRAS  2055 1-4250 .  4.0  x  10“6  3.5  x  10“4  0.011 

IRAS  23365+3604 .  1.5  x  10“6  2  x  10~4  0.0075 


Normals/Starbursts 

M82 .  1 .1  x  10-3  a  4  x  10  1  2.8 

CenA .  3.3  x  10_4b  2.5  x  10“4  1.3 

NGC  1068 .  2.5  x  10_4c  2.5  \  10  4  1.0 

NGC1313 .  1.4  x  10_5d  3.5  x  II)  4  0.04 

NGC3256 .  3.8  x  10_4e  3.0xl0~4  1.3 

NGC  3690  (Arp  299) .  1.0  xl0-4f  3.0xl0“4  0.33 


a  Colbert  et  al.  1999. 
b  Unger  et  al.  2000. 
c  L.  Spinoglio  et  al.  (in  preparation). 
d  Contursi  et  al.  2002  (region  90). 
e  Carralet  al.  1994. 
f  Satyapal  et  al.  2003. 


predictions  and  previous  observations  of  CO  (1-0),  we  can 
estimate  the  CO  ( 1-0)  beam  area  filling  factor  $co  as  we  did 
with  [C  ii]. 

Table  4  gives  the  observed  CO  (1-0)  intensities  for  Mrk 
273,  IRAS  17208-0014,  IRAS  20551-2450,  and  IRAS 
23365+3604  from  previous  studies  scaled  to  the  69"  beam 
of  our  ISO  [C  ii]  data.  Again,  taking  the  predicted  n  and  G0 


for  these  galaxies  from  §  3.1,  we  list  in  Table  4  the  predicted 
CO  (1-0)  intensities  according  to  Kaufman  et  al.  (1999). 
The  ratio  of  the  measured  to  predicted  CO  (1-0)  intensities 
yields  the  CO  (1-0)  beam  area  filling  factor  <I\  ().  Comparing 
the  CO  and  [C  ii]  beam  area  filling  factors,  we  find  that  $co 
is  4  to  1 1  times  larger  than  $[C  „]  (Table  4,  col.  [5]).  In  other 
words,  within  the  ISO  69"  beam,  the  effective  surface  area  of 


TABLE  4 

CO  Beam  Area  Filling  Factors  <bco  for  Four  ULIRGs  and  Comparison  Normal/Starburst  Galaxies 


Observed  Go  (i-0)a  Predicted  Go  (1-0)  f°r  <&co  =  1 
Source  (ergscm_2s_1sr_1)  (ergscm_2s_1sr_l)  4>co  'f’co/'t’p  n] 

(1)  (2)  (3)  (4)  (5) 


ULIRGs 

Mrk  273 .  2.4  xl0-9b  2  x  10“8  0.12  4.0 

IRAS  17208-0014 .  4.6  xl0-9c  2  x  10“8  0.23  5.9 

IRAS 2055 1-4250 .  2.0  xl0-9d  4.5  x  10'8  0.044  4.0 

IRAS  23365+3604 .  1.7  xl0-9d  2  x  10“8  0.085  11 


Normals/Starbursts 

M82 .  i .7  x  HI  7e  6.5  x  10  8  2.6  0.9 

CenA .  7.6  x  10"8 e  5.0  x  10'8  1.5  1.2 

NGC  1068 .  5.8  xl0-8e  5.3  x  10~8  1.1  1.1 

NGC  1313 .  4.1  x  10_9f  7.5  xlO-8  0.05  1.3 

NGC  3256 .  1.0  x  10~7g  5.5  x  10'8  1.8  1.4 

Arp  299 .  1.5  x  10~8  h  4.0  x  10'8  0.37  1.1 


a  The  CO  intensities  are  scaled  to  the  ISO  [C  ii]  beam  size  (69").  For  galaxies  that  are  small  compared  to  the 
beam  size,  we  divide  the  observed  CO  intensity  by  8.8  x  10-8  sr  (69");  for  galaxies  larger  than  the  beam  size,  we 
multiply  the  observed  CO  intensity  by  the  ratio  of  the  ISO  beam  size  to  the  beam  size  of  the  CO  observations. 
b  Sanders,  Scoville,  &  Soifer  1991  (55"  beam). 
c  Mirabel  et  al.  1990  (44"  beam). 
d  Solomon  et  al.  1997  (22"  beam). 
e  Stacey  et  al.  1991  (55"  beam). 
f  Contorsi  et  al.  2002  (43"  beam), 
s  Aalto  et  al.  1995  (43"  beam). 
h  Sanders  &  Mirabel  1985  (60"  beam). 
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the  environment  that  gives  rise  to  CO  ( 1-0)  emission  is  4  to 
1 1  times  that  of  the  gas  that  emits  [C  n]. 

Having  established  that  the  [C  ii]  and  CO  filling  factors 
differ  significantly  for  the  four  ULIRGs  in  Tables  3  and  4, 
and  possibly  for  ULIRGs  in  general,  the  question  now 
becomes:  How  does  4>[C  n]  compare  to  the  FIR  beam  area 
filling  factor  $FIR?  We  cannot  use  the  approach  above  to 
independently  derive  <bFIR  for  comparison  with  $[C  „],  since 
both  the  observed  [C  n]  and  FIR  flux  were  used  to  derive  the 
PDR  n  and  G0.  Instead,  to  address  this  question,  we  refer  to 
a  review  paper  by  Genzel  &  Cesarsky  (2000).  In  Figure  8 a  of 
their  paper,  Genzel  &  Cesarsky  plot  the  intensity  ratio 
%:ii]/4ir  against  Ico  (i_o)//Fir  for  a  variety  of  galaxy  types 
including  ULIRGs.  Using  PDR  model  overlays  (Wolfire, 
Hollenbach,  &  Tielens  1989)  to  provide  a  measure  of  Go  and 
n,  Genzel  &  Cesarsky  find  that  the  location  of  the  ULIRGs 
in  their  plot  immediately  suggests  that  the  [C  ii]  deficit  can 
be  quantitatively  understood  if  the  gas  densities  exceed  105 
cm~3,  with  Go  ~  103  104,  typical  of  Galactic  PDR/H  ii  star 
forming  regions.  This  finding  assumes,  however,  that  the 
[C  ii],  CO,  and  FIR  beam  area  filling  factors  are  identical.  In 
this  proposed  high-density  case  where  n  lies  well  above  the 
[C  ii]  critical  density  (3  x  103  cm~3  for  collisions  with  hydro¬ 
gen;  Flower  &  Launay  1977;  Hayes  &  Nussbaumer  1984), 
the  abundance  of  C+  in  the  upper  level  is  fixed  by  the  atomic 
statistical  weights.  Furthermore,  the  column  density  of  C+ 
becomes  insensitive  to  G0  (Wolfire,  Hollenbach,  &  Tielens 
1989).  Since  the  abundance  and  column  density  cease  to  rise 
with  n  and  G0,  the  emergent  [C  ii]  line  intensity  saturates, 
giving  rise  to  a  [C  ii]  deficit  with  respect  to  the  FIR  emission. 

The  interpretation  that  the  [C  ii]  line  arises  in  high-density 
PDR/H  ii  regions  cannot  adequately  explain  the  PDR  line 
emission  for  ULIRGs,  however.  In  particular,  according  to 
Figure  6,  we  see  that  for  densities  of  105  cm~3  and 
Go  ~  103  104,  the  strength  of  the  [O  i]  145  /in l  line  should  be 
comparable  to  that  of  the  [C  ii]  line.  The  observed  [O  i]  145 
/an  line  strengths  in  the  ULIRGs,  on  the  other  hand,  fall 
well  below  that  of  the  observed  [C  ii]  line.  In  addition, 
assuming  no  [O  i]  63  /an  self-absorption,  the  Kaufman  et  al. 
PDR  models  imply  a  [O  i]  63  /an  to  [C  n]  ratio  of  10-30  for 
densities  of  105  cnr3  and  Go  ~  103  104,  which  is  much 
higher  than  observed.  This  mismatch  with  the  [O  i]  observa¬ 
tions  implies  again  that  PDR  models,  together  with  the 
assumption  that  [C  ii],  CO  (1-0),  and  FIR  beam  filling 
factors  are  similar,  are  incompatible  with  the  observations. 

We  have  shown  that  for  the  four  ULIRGs  for  which  we 
have  both  [C  ii]  and  unabsorbed  [O  i]  63  /an  emission  line 
fluxes,  4>co  is  4  to  11  times  larger  than  <b[C  „].  If  we  assume 
that  <h[ClI]  =  4>fir,  i.e.,  that  the  observed  [C  n]  to  FIR  ratio 
is  the  intrinsic  ratio  unaffected  by  beam  area  filling  factor 
effects,  then  the  observed  ratio  of  /Co  (i-o)Afir  must  be  div¬ 
ided  by  a  factor  of  4  to  11  (=4)co/‘j)[c  n]  or  'bco/  $fir) 
before  plotting  the  ULIRG  data  in  Figure  8a  (/[Ch]//fir  vs. 
/Co  (i-o)/^fir)  of  Genzel  &  Cesarsky  (2000).  We  find  that 
such  a  correction,  however,  would  result  in  only  a  factor  of 
a  few  change  in  the  derived  Go  and  n,  and  these  values  would 
still  remain  inconsistent  with  the  observed  [C  ii]  and  [O  i] 
63  /mi  and  145  //m  line  ratios. 

On  the  other  hand,  if  4>Co  =  ^fir  for  the  ULIRGs,  then 
4>[ClI]  and  4>FiR  must  differ,  and  the  observed  /[Cii]/^fir  is 
not  the  intrinsic  ratio.  Specifically,  for  the  four  galaxies  in 
Tables  3  and  4,  we  would  correct  /[oi]/7fir  by  multiplying 
by  <f>co/<b[c  ii]  (=$fir/$[c  h]  for  ‘bco  =  $fir)  in  Table  3, 
column  (5).  In  this  instance,  the  location  of  the  ULIRGs  in 


Figure  8a  of  Genzel  &  Cesarsky  implies  a  gas  density 
n  ~  103  104  cm~3  and  a  UV  flux  Go  ~  103  104.  Such 
derived  PDR  properties  are  typical  of  normal  and  starburst 
galaxies  and  would  provide  agreement  between  the 
observed  [O  i]  63  /in i  and  145  //m  and  [C  ii]  line  ratios  and 
the  PDR  models.  Furthermore,  such  a  beam  area  filling 
factor  correction  would  remove  the  [C  ii]  deficit,  yielding  a 
[C  ii]  to  FIR  ratio  closer  to  0.1%.  Thus,  we  could  account  for 
the  [C  ii]  deficit  in  ULIRGs  with  Go  and  n  parameters  typical 
for  starburst  galaxies  if  the  beam  area  filling  factor  for  FIR  is 
significantly  larger  than  that  for  [C  ii],  as  it  is  for  CO.  Here, 
the  PDR  model  comparisons  that  support  higli-Go/n  PDRs  in 
ULIRGs  would  not  apply  because  the  beam  area  filling  factors 
for  [C  ii]  ( and  presumably  [O  i])  and  FIR  are  not  the 
same.  The  presence  of  additional  non-PDR  components 
of  FIR  and  CO  emission  would  provide  one  possible 
means  of  achieving  FIR  and  CO  luminosity  without  the 
accompanying  [C  ii]  emission. 

If  the  above  analysis  is  correct,  we  should  see  agreement 
in  the  beam  area  filling  factors  for  [C  ii]  and  CO  in  normal 
and  starburst  galaxies,  in  contrast  to  the  ULIRGs.  To  test 
this  condition,  we  performed  the  above  beam  area  filling 
factor  analysis  for  an  assortment  of  normal  and  starburst 
galaxies  for  which  there  exist  [C  ii]  and  [O  i]  data  (using  ISO 
data  where  available)  and  CO  (1-0)  data  with  spatial  resolu¬ 
tion  comparable  to  the  [C  ii]  and  [O  i]  observations.  The 
[C  ii]  and  CO  beam  area  filling  factors  for  these  galaxies  are 
shown  in  Tables  3  and  4.  In  all  cases,  3>co/®[Cn]  is  near 
unity  for  the  normal  and  starburst  galaxies,  in  contrast  to 
our  finding  for  the  ULIRGs  (Table  4,  col.  [5]).  [We  note  that 
for  some  normal  galaxies,  particularly  early-type  galaxies, 
characterized  by  low  star  formation  rates,  the  [C  ii]  emission 
may  arise  primarily  from  the  diffuse  ISM,  in  which  case  we 
would  not  expect  4>[C  „  and  <I>(  0  to  agree  (Pierini  et  al. 
2001).]  Thus,  we  conclude  that,  for  ULIRGs,  the  beam  area 
filling  factor  differences  between  [C  ii]  and  CO  are  real. 

The  presence  of  CO  emission  arising  largely  from  a  non- 
PDR  component  where  there  is  no  [C  ii]  is  consistent  with 
CO  observations  of  Arp  220  and  Mrk  231  (e.g.,  Downes  & 
Solomon  1998;  Scoville  et  al.  1997).  These  observations 
show  that  much  of  the  CO  emission  arises  from  molecular 
gas  that  is  not  cospatial  with  the  nuclei  (the  presumed  sites 
of  the  FIR  emission)  as  seen  in  the  millimeter  continuum. 
The  suggestion  that,  like  CO  (1-0),  the  FIR  emission  arises 
in  large  part  from  a  non-PDR  component  provides  a  com¬ 
pelling  explanation  for  the  [C  ii]  deficit,  although  without 
the  benefit  of  an  FIR  interferometer,  conclusive  evidence 
that  this  is  the  case  is  difficult  to  produce.  Although  the  loca¬ 
tion  of  the  [C  ii]  emission  is  uncertain,  we  note  that  it  is  not 
unreasonable  to  suggest  that  the  [C  ii]  emission  may  be  coin¬ 
cident  with  the  PAH  emission  (discussed  further  in  §  3.2.3), 
which  Soifer  et  al.  (2002)  have  found  is  somewhat  extended 
and  not  fully  coincident  with  the  mid-IR  continuum  in 
Arp  220  and  Mrk  273.  This  would  tend  to  support  the  idea 
that  the  FIR  and  [C  ii]  emission  may  not  be  cospatial. 

3.2.2.  C  i  (1—0)  Emission  Component 

The  609  /im  [C  i]  (1-0)  line  is  another  fine-structure  line 
arising  in  PDRs  for  which  there  exist  ULIRG  observations. 
Again,  these  observations  and  the  PDR  models  allow  us  to 
derive  a  [C  i]  beam  area  filling  factor  $[C  i]  for  comparison 
with  <I» rC  nj .  According  to  one-dimensional  PDR  models 
(e.g.,  Tielens  &  Hollenbach  1985),  in  clouds  exposed  to  UV 
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radiation,  carbon  is  mostly  in  ionized  form  to  a  depth 
Ay  ~  1  mag.  Atomic  carbon  appears  at  an  intermediate 
(A  y  ~  1-5  mag)  depth,  where  C+  has  mostly  recombined 
with  electrons  to  form  C  and  not  all  of  the  gas-phase  carbon 
is  converted  into  CO.  The  [C  i]  line  thus  traces  regions 
exposed  to  nonionizing  UV  radiation.  Since  FIR  continuum 
radiation  comes  from  regions  exposed  to  both  ionizing  and 
nonionizing  UV  radiation,  estimates  of  $[C  ,]  can  place  a 
lower  limit  on  4>fir. 

Gerin  &  Phillips  (2000)  measure  the  [C  i]  (1-0)  emission 
in  a  sample  of  12  galaxies,  in  which  two  galaxies  are 
ULIRGs — Arp  220  and  Mrk  23 1 .  In  the  case  of  Arp  220,  we 
have  good  estimates  of  n  and  G0  from  previous  studies, 
which  we  can  use  to  determine  4>[(- Gerin  &  Phillips  meas¬ 
ure  a  [C  i]  (1-0)  intensity  of  2.16  x  1CU7  ergs  cm~2  s~!  sr~* 
scaled  to  the  ISO  69"  beam.  CO  observations  of  Arp  220 
suggest  a  moderately  dense  interstellar  medium  with 
n  ~  104  cm  3  (Solomon  et  al.  1997;  Scoville  et  al.  1997; 
Gerin  &  Phillips  1998).  Based  on  the  FIR  continuum  and 
hydrogen  recombination  line  emission  in  Arp  220,  Gerin  & 
Phillips  (1998)  estimate  the  UV  radiation  field  Go  to  be 
4  x  104.  According  to  Kaufman  et  al.  (1999),  the  predicted 
[C  i]  (1-0)  intensity  assuming  unity  beam  area  filling  factor 
is  5.1  x  10-6  ergs  cm~2  s_1  sr_1.  The  ratio  of  the  observed 
[C  i]  intensity  to  the  predicted  intensity  yields  a  [C  i]  (1-0) 
beam  area  filling  factor  4>[C ,j  of  0.042  for  a  69"  beam.  Simi¬ 
larly,  taking  the  [C  n]  intensity  in  Arp  220  (9.9  x  10~6  ergs 
cm~2  s  1  sr  ! )  and  comparing  with  the  [C  ii]  intensity  pre¬ 
dictions  of  Kaufman  et  al.  (1999)  for  a  unity  beam  area  fill¬ 
ing  factor,  we  derive  a  [C  ii]  beam  area  filling  factor  <I>|C  nj  of 
0.01  for  Arp  220,  or  <f?[Cl]/‘£[cII]  ~  4.  For  comparison,  using 
the  unity  beam  area  filling  factor  CO  intensity  as  predicted 
by  Kaufman  et  al.,  we  estimate  that  <bco  =  0.17 
Gcofi-o)  =  2.05  x  1(U8  ergs  cm~2  s_1  sr~'  corrected  to  a  69" 
beam;  Radford,  Solomon,  &  Downes  1991;  Gerin  &  Phillips 
1998),  or  $co/4>[ClI]  ~  17. 

Thus,  as  with  the  CO  (1-0)  emission  in  Mrk  273,  IRAS 
17208-0014,  IRAS  20551-4250,  and  IRAS  23365+3604, 
we  see  that  the  implied  projected  area  of  the  [C  i]-emitting 
region  (and  CO-emitting  region)  is  significantly  larger  than 
that  of  the  [C  n]-emitting  region  in  Arp  220.  We  contend,  as 
stated  above,  that  the  value  of  4>[C  ,]  places  a  crude  lower 
limit  on  <1>F|R.  Multiplying  the  observed  value  of  F]Cn]/.FFIR 
in  Arp  220  by  $Fir/ 4>|C  uj  can  therefore  account  for  most,  if 
not  all,  of  the  [C  ii]  deficit. 

For  the  sake  of  argument,  let  us  assume  that  the  [C  ii] 
emission  from  PDRs  is  coextensive  with  the  [C  i]  and  FIR 
emission,  i.e.,  that  the  [C  ii]  is  intrinsically  depressed  with 
respect  to  [C  i]  and  FIR.  Discounting  opacity  and  extinction 
effects  (discussed  in  §  3.4),  PDR  models  should  therefore 
apply,  and  a  unique  G0  and  n  regime  that  explains  the  [C  n], 
[C  i],  and  FIR  emission  should  exist.  For  the  12  galaxies  in 
their  sample,  Gerin  &  Phillips  (2000)  plot  [C  n]/FIR  and 
[C  ii]/[C  i]  versus  [C  i]/FIR  on  contours  of  n  and  G0  from 
PDR  model  predictions.  Most  of  the  galaxies  studied  lie  in 
the  parameter  space  Go  =  102— 103,  n  =  102  105,  typical  of 
molecular  clouds  in  the  nuclei  of  nearby  galaxies.  However, 
Arp  220  and  Mrk  231 — the  only  ULIRGs  in  the  sample — 
lie  outside  of  the  parameter  space  depicted  by  the  modeled 
contours,  but  in  the  direction  of  very  low  G0/n  (s£0.01).  Such 
values  of  G0/n  are  at  least  2  orders  of  magnitude  lower  than 
the  estimates  obtained  with  the  FIR  fine-structure  lines, 
again  suggestive  of  filling  factor  discrepancies  attributable 
to  a  non-PDR  component(s). 


3.2.3.  6.2  //m  PAH  Emission  Component 

Referring  again  to  Figure  2,  the  depicted  [C  n]  to  PAH 
emission  behavior  provides  compelling  support  for  the 
influence  of  filling  factor  effects  on  the  observed  emission.  If 
the  above  filling  factor  analysis  is  correct,  we  should  measure 
normal  line  ratios  when  comparing  [C  n]  with  another  species 
that  is  known  to  come  from  the  same  gas.  Fike  [C  n],  the 
6.2  //m  PAH  emission  traces  gas  exposed  to  far-UV  radiation 
(Tielens  &  Hollenbach  1999;  Vigroux  et  al.  1999).  In  Figure 
2,  we  see  that  in  fact  the  relative  [C  ii]  to  PAH  emission 
remains  unchanged  when  comparing  the  ULIRGs  to  the 
normal  and  starburst  galaxies.  We  would  expect  the  observed 
difference  in  the  [C  ii]  to  FIR  and  6.2  /on  PAH  to  FIR  ratios 
between  ULIRGs  and  normal/starburst  galaxies  to  be  a 
normal  by-product  of  an  additional  FIR  component  in 
ULIRGs  that  is  not  found  in  less  luminous  galaxies.  Soifer 
et  al.  (2002)  have  confirmed  that  the  mid-IR  continuum  in 
Arp  220  and  Mrk  273  is  strongly  absorbed,  while  showing  for 
the  first  time  that  the  11.3  /mi  PAH  emission  is  extended.  If 
the  absorbed  mid-IR  radiation  is  powering  the  FIR  emission, 
then  even  an  apparent  overlap  between  the  PAH  emission 
and  the  FIR  emission  would  not  mean  that  they  share  the 
same  power  source. Thus,  within  the  ISO  LWS  beam,  as  with 
CO  and  [C  i],  there  could  likely  exist  sources  of  FIR  emission 
lacking  associated  PAH  and  [C  ii]  emission  that,  once  taken 
into  account,  could  explain  the  [C  ii]  deficit. 

3.2.4.  Implications  for  Derived  PDR  Parameters 

Although  we  are  unable  to  derive  precise  corrected  PDR 
parameters  resulting  from  the  inferred  presence  of  a  non- 
PDR  FIR  component,  we  note  that  such  a  correction  would 
move  the  ULIRGs  to  the  right  in  Figure  4,  resulting  in 
higher  densities  and  lower  Go  values  compared  with  the 
uncorrected  values.  Depending  on  the  amount  of  the  correc¬ 
tion,  this  is  more  consistent  with  the  higher  densities 
thought  to  be  present  in  ULIRGs  (Solomon  et  al.  1997).  In 
this  case,  the  warm  60-100  //m  colors  of  ULIRGs  would  be 
generated  in  the  non-PDR  component  of  the  FIR  emission. 
We  can  see  from  Figure  6  that  parameter  changes  in  these 
directions  (lower  Go  and  higher  n)  would  not  have  a  strong 
effect  on  the  [C  ii]  to  [O  i]  145  /mi  ratio. 

3.3.  Soft  UV  Radiation  Field:  Aging  Starbursts 

In  Paper  I,  a  soft  UV  radiation  field  was  presented  as  a 
possible  explanation  for  the  [C  ii]  deficit  in  the  observed 
ULIRGs.  A  stellar  population  deficient  in  massive  main- 
sequence  stars,  as  a  result  of  either  an  aging  starburst  or  an 
initial  mass  function  with  a  low  upper  mass  cutoff,  could 
give  rise  to  a  soft  UV  field,  and  consequently  a  small  Fy  n]  / 
Frr  ratio  (Spaans  et  al.  1994;  Nakagawa  et  al.  1995; 
Malhotra  et  al.  2001).  S.  Satyapal  et  al.  (in  preparation) 
model  the  decrease  in  L[Ch]/Lfir  as  a  function  of  starburst 
age.  In  their  model,  an  ionization-bounded  H  ii  region  and 
PDR  interface  surround  a  central  starburst  cluster,  with  the 
location  of  the  PDR  gas  lying  at  the  edge  of  the  ionized 
region.  As  the  radiation  field  softens,  there  is  a  reduction  in 
the  UV  photons  that  heat  the  gas,  resulting  in  a  decrease  in 
the  [C  ii]  flux.  Since  both  UV  and  visible  photons  can  heat 
the  grains,  the  [C  ii]  to  FIR  ratio  decreases  with  age. 

We  do  not  favor  a  soft  UV  field  arising  from  an  aging 
starburst  as  the  explanation  for  the  [C  ii]  deficit.  Since  star- 
burst  luminosities  decrease  as  the  stellar  population  ages, 
older  starbursts  in  ULIRGs  would  imply  that  these  same 
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ULIRGs  were  much  more  luminous  in  the  past,  inconsistent 
with  the  observed  luminosity  distribution  of  1R  galaxies  in 
the  local  universe.  This  explanation  is  therefore  not  very 
compelling,  since  it  would  imply  that  the  less  luminous 
objects  are  characterized  by  younger  starbursts.  S.  Satyapal 
et  al.  (in  preparation)  model  the  [C  n]  to  FIR  ratio  as  a  func¬ 
tion  of  starburst  age.  They  show  that  if  an  aging  starburst 
produces  the  low  [C  n]  to  FIR  ratio  in  Arp  220,  the  starburst 
taking  place  in  this  galaxy  would  need  to  be  at  least  3  x  106 
yr  older  than  the  one  taking  place  in  the  less  luminous  gal¬ 
axy  Arp  299.  If  the  starburst  in  Arp  220  were  3  x  106  yr 
older  than  the  one  taking  place  in  Arp  299,  the  luminosity  of 
Arp  220  would  have  been  roughly  3  times  greater  than  its 
current  value  several  million  years  ago  and  almost  an  order 
of  magnitude  greater  a  couple  of  million  years  earlier.  If 
UFIRGs  are  indeed  older  starbursts,  we  would  expect  to  see 
a  population  of  >1013  L0  UFIRGs  in  the  local  universe, 
which  we  do  not. 

3.4.  Optical  Depth  and  Extinction  Effects 

In  Paper  I,  we  mentioned  the  possibility  of  self-absorbed 
or  optically  thick  [C  ii]  emission.  As  discussed,  self-absorbed 
[C  ii]  requires  that  C+  exist  in  cooler  foreground  material. 
Self-absorbed  [C  n]  emission  has  been  observed  toward  the 
H  ii  region/molecular  cloud  complex  Sgr  B2  near  the  Galac¬ 
tic  center  (Coxet  al.  1999;  Vastel  et  al.  2002).  In  Sgr  B2,  FIR 
fine-structure  line  emission  from  ionized  and  PDR  gas  is 
extremely  faint  (Cox  et  al.),  as  is  the  case  in  UFIRGs,  and 
the  [O  i]  63  pm  line  is  seen  in  absorption,  as  is  the  case  in 
some  UFIRGs,  such  as  Arp  220.  However,  we  do  not 
believe  that  the  [C  n]  deficit  observed  in  UFIRGs  is  caused 
primarily  by  self-absorption.  Rather,  we  suggest  (and  dis¬ 
cuss  further  in  the  next  section)  that  the  cause  of  the  [C  ii] 
deficit  in  these  sources  is  also  responsible  for  the  faintness  of 
the  [O  i]  145  pm  line,  which  we  note  does  not  appear  to  show 
self-absorption  in  Sgr  B2  (Vastel  et  al.  2002),  as  well  as  that 
of  the  [O  iii]  52  and  88  pm  and  [N  hi]  57  pm  lines.  As  for 
optically  thick  [C  n],  the  C+  column  densities  needed  to 
attain  the  required  optical  depth  in  the  [C  ii]  line,  i.e., 
-7-[ClI]  5=2,  correspond  to  an  unrealistically  high  Ay  (>130 
mag)  in  the  PDR  gas  alone  compared  to  previous  measures 
of  the  total  Av  for  the  cold  molecular  and  PDR  gas  (Sturm 
et  al.  1996). 

Comparison  of  mid-IR  recombination  and  fine-structure 
line  fluxes  with  the  free-free  millimeter  continuum  in  the 
prototypical  UFIRG  Arp  220  suggests  that  extinction  high 
enough  to  produce  obscuration  at  FIR  wavelengths,  or 
specifically  158  pm ,  is  not  evident  (Genzel  et  al.  1998; 
S.  Satyapal  et  al.,  in  preparation),  at  least  for  the  case  in 
which  the  gas  responsible  for  the  observed  mid-IR  recombi¬ 
nation  and  fine-structure  line  emission  is  also  responsible 
for  the  millimeter  emission.  This  result  would  tend  to  rule 
out  high  extinction  at  158  pm  as  the  origin  of  the  [C  ii]  158 
pm  deficit.  Furthermore,  significant  extinction  at  158  pm  is 
not  supported  by  the  relatively  normal  [O  i]  63  pm  to  [C  ii] 
158  //nr  line  ratios  in  the  UFIRGs  in  this  study  and  the  ten¬ 
dency  toward  higher  [O  i]  to  [C  ii]  ratios  in  the  [C  ii]  deficient 
galaxies  in  the  study  of  Malhotra  et  al.  (2001),  since  high  dif¬ 
ferential  extinction  would  act  to  decrease  the  [O  i]  63  pm 
emission  relative  to  [C  n]  and  produce  less  than  normal  [O  i] 
63  pm  to  [C  n]  line  ratios.  For  A\ooA^a  and  l>a=F2,  we 
would  need  to  multiply  the  observed  [O  i]  63  //nr  to  [C  ii]  line 
ratio  by  a  factor  of  at  least  16  (for  a=  1,  and  at  most 


2.5  x  105  for  a  =  2)  to  correct  for  extinction  if  extinction 
at  158  pm  alone  accounts  for  the  [C  n]  deficit,  i.e.,  if 
^158/im  ~  2.  Such  a  correction  to  the  observed,  relatively 
normal  [O  i]  to  [C  ii]  ratios  would  generate  flux  ratios  greatly 
in  excess  of  the  predictions  of  the  PDR  models  for  any 
reasonable  n  and  Go- 

Although  extinction  at  FIR  wavelengths  cannot  explain 
the  [C  ii]  deficit,  mid-IR  spectroscopic  observations  by 
Genzel  et  al.  (1998)  do  appear  to  indicate  higher  nrid-IR 
extinction  toward  the  ionized  medium  in  UFIRGs  than  in 
normal  and  starburst  galaxies.  Because  of  uncertainties  in 
the  nature  and  exact  shape  of  the  extinction  curve  in  the 
mid-IR,  the  degree  to  which  extinction  may  lower  the  PAH 
flux  relative  to  FIR  in  UFIRGs  as  compared  to  normal  and 
starburst  galaxies  is  difficult  to  quantify.  For  example,  the 
Draine  (1989)  mid-IR  extinction  law  (A^ooA^1  -75)  implies  a 
factor  of  only  1.7  drop  in  the  6.2  pm  PAH  flux  going  from 
Ay  =  10  (M82;  Genzel  et  al.  1998)  to  Av  =  45  (Arp  220; 
Genzel  et  al.  1998).  To  explain  the  PAH  deficit  relative  to 
FIR  in  UFIRGs  would  require  Ay  >  100,  in  disagreement 
with  the  extinction  estimates  of  mid-IR  spectroscopy.  Futz 
(1999),  however,  derives  a  mid-IR  extinction  curve  based  on 
ISO  observations  of  the  highly  obscured  Galactic  center 
region  that  indicates  significantly  higher  extinction  in  the  3- 
8  pm  region  than  implied  by  Draine.  According  to  Futz,  the 
6.2  pm  flux  should  drop  by  nearly  a  factor  of  5,  going  from 
Av  =  10  to  45,  in  which  case  mid-IR  extinction  could 
account  for  much  of  the  PAH  deficit  if  the  PAH  emitting 
regions  suffer  the  same  extinction  as  that  derived  from  the 
recombination  and  fine-structure  lines.  However,  this  may 
not  be  the  case.  For  example,  recent  evidence  suggests  that 
the  mid-IR  continuum  in  AGNs  may  suffer  more  extinction 
than  the  PAHs.  Clavel  et  al.  (2000)  find  that,  although  the 
mean  equivalent  width  of  the  7.7  pm  feature  in  Seyfert  1  gal¬ 
axies  is  a  factor  of  5.4  lower  than  that  of  their  Seyfert  2 
counterparts,  the  mean  PAH  luminosities  of  the  two  types 
are  equal.  Clavel  et  al.  attribute  the  high  mean  equivalent 
width  in  Seyfert  2  galaxies  to  extinction  of  the  mid-IR  con¬ 
tinuum.  This  would  imply  that  in  Seyfert  1  and  2  galaxies, 
the  PAH  extinction  is  the  same  even  though  the  continuum 
emission  suffers  very  different  extinction. 

Given  the  uncertainties  in  the  mid-IR  extinction,  the  rela¬ 
tive  constancy  in  the  PAH  to  [C  ii]  ratio  between  UFIRGs 
and  normal/starburst  galaxies,  as  well  as  in  the  Helou  et  al. 
(2001)  sample,  is  indeed  remarkable.  It  would  therefore  be 
somewhat  surprising  if  effects  such  as  lower  photoelectric 
heating  (which  would  lower  the  [C  n]  to  FIR  flux)  and 
higher  average  differential  extinction  (which  would  lower 
the  PAH  to  FIR  flux)  in  UFIRGs  just  cancel  each  other  to 
produce  the  observed  constant  ratio,  especially  considering 
the  large  “  moment  arm  ”  in  wavelength  space  between  the 
6.2  //nr  PAH  feature  and  the  [C  ii]  158  pm  line. 

3.5.  Dust  within  Photoionization  Environments 

Based  on  the  analysis  of  §  3.2,  we  see  that  the  application 
of  PDR  models  under  simplifying  assumptions  does  not 
yield  consistent  results.  The  differences  in  [C  ii],  CO,  [C  i], 
and  presumably,  FIR  beam  area  filling  factors  imply  that  a 
significant  fraction  of  the  FIR  could  come  from  environ¬ 
ments  where  there  is  little  or  no  [C  ii]  or  PAH  emission. 
Photoionization  regions  may  provide  one  such  environment 
if  the  dust  in  the  ionized  gas  is  a  primary  contributor  to  the 
overall  FIR  luminosity. 
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In  photoionization  regions,  dust  grains  compete  with  the 
gas  in  absorbing  the  direct  stellar  and  nebular  continuum 
radiation.  Dust  grains  reradiate  the  absorbed  photons  as 
thermal  continuum  radiation  in  the  IR.  We  would  expect 
ionized  gas  with  dust  grains  that  compete  efficiently  with  the 
hydrogen  for  the  available  UV  photons  to  exhibit  enhanced 
FIR  emission  relative  to  the  recombination  and  fine- 
structure  line  emission  front  the  gas.  Compared  to  their  less 
luminous  counterparts,  the  mid-  to  far-IR  spectra  of 
ULIRGs  do  in  fact  display  a  lack  of  strong  fine-structure 
emission  lines  from  ionized  gas,  as  well  as  an  abundance  of 
molecular  absorption  lines  (Lutz  et  al.  1996;  Genzel  et  al. 
1998;  Fischer  et  al.  1999).  Fischer  et  al.  argue  that  the  weak¬ 
ness  of  the  FIR  fine-structure  lines  is  not  due  to  FIR  extinc¬ 
tion.  Voit  (1992)  suggests  that  direct  absorption  of  the 
ionizing  luminosity  by  dust  could  reduce  fine-structure  and 
recombination  line  emission  in  ULIRGs.  The  warmer  FIR 
colors  and  dust  temperatures  of  ULIRGs  (e.g.,  Sanders  et 
al.  1988)  compared  to  less  luminous  galaxies  suggest  that 
the  role  of  dust  in  ionized  gas  is  important. 

As  discussed  by  Voit  (1992)  and  Bottorff  et  al.  (1998), 
dust  absorbs  most  of  the  ionizing  photons  when  the  dimen¬ 
sionless  ionization  parameter  U  is  high,  where  U  is  defined 
as 

jj  _  Q  _  $(H)  _  Niinas 

4ir  r2nec  nec  c 

where  Q  and  $(11)  are  the  rate  and  flux  of  ionizing  photons, 
respectively,  ne  is  the  electron  density,  c  is  the  speed  of  light, 
Ah  ii  is  the  ionized  hydrogen  column  density,  and  aB  is  the 
case  B  hydrogen  recombination  coefficient  (Storey  & 
Hummer  1995).  Extragalactic  clouds  with  larger  ionization 
parameters  have  larger  ionized  column  densities  and  larger 
grain  column  densities,  which  results  in  the  absorption  by 
dust  of  a  larger  fraction  of  the  incident  luminosity.  For 
clouds  with  standard  Galactic  dust  abundance  and  compo¬ 
sition,  dust  absorbs  most  of  the  incident  ionizing  radiation 
above  U  =  1(U2  (Voit  1992;  Bottorff  et  al.  1998).  In  such 
cases,  the  ionized  hydrogen  column  density  is  typically  1021 
cm~2  or  greater.  Bottorff  et  al.  show  that  the  IR  to  H (3  ratio 
is  a  strong  function  of  U,  with  the  IR  intensity  increasing  as 
grains  absorb  more  radiation  and  H/3  decreasing  as  hydro¬ 
gen  absorbs  fewer  ionizing  photons.  Such  “  dust-bounded  ” 
cases  fall  in  the  IR/H/3^>100  limit  (Bottorff  et  al.  1998). 
Based  on  models  with  the  photoionization  code  CLOUDY 
(Ferland  1993),  Fischer  et  al.  (2001)  find  that  high-U condi¬ 
tions  may  explain  the  lack  of  detected  fine-structure  line 
emission  in  ULIRGs.  For  example,  these  models  predict 
that  for  central  instantaneous  starbursts  and  power-law 
sources  with  Q  =  4.5  x  1054  s_1  and  ne  =  500  cm~3, 
lowering  the  distance  between  the  central  source  and  the 
surrounding  gas  from  500  pc  (U  =  10~2)  to  50  pc  (U  =  10°) 
decreases  the  ratio  of  the  [O  iii]  88  /mi  line  to  the  total  lumi¬ 
nosity  by  more  than  an  order  of  magnitude. 

Enhanced  dust  absorption  in  the  ionized  gas  offers  a 
natural  explanation  for  the  implied  beam  area  filling  factor 
differences  between  the  [C  n]  emission  and  the  CO, 
[C  i],  and  FIR  emission.  If  a  significant  portion  of  the  FIR 
emission  emanates  from  the  photoionized  gas,  while  the 
[C  n]  does  not,  we  would  expect  this  extra  source  of  FIR 
emission  to  increase  the  apparent  FIR  beam  area  filling  fac¬ 
tor  relative  to  that  of  [C  n].  Here  again,  we  emphasize  that 
unless  the  FIR  emission  from  the  ionized  gas  is  taken  into 


account,  either  by  applying  a  beam  area  filling  factor  correc¬ 
tion  or  by  subtracting  the  FIR  contribution  from  the  ion¬ 
ized  gas  from  the  total  FIR  emission,  the  PDR  model 
predictions  in  §  3.1  do  not  apply  to  ULIRGs.  By  inhibiting 
the  penetration  of  13.6-6  eV  photons  (i.e.,  those  photons 
thought  to  be  responsible  for  heating  the  PDR  gas)  and 
decreasing  the  extent  of  the  [C  ii]-emitting  region,  dust- 
bounded  ionized  gas  should  also  give  rise  to  [C  ii]  beam  area 
filling  factors  that  are  smaller  than  those  for  CO  and  [C  i]. 
This  is  indirectly  supported  by  numerous  extragalactic 
observations  of  an  effect  in  the  opposite  sense,  whereby  a 
reduced  dust  abundance  allows  the  UV  radiation  to  pene¬ 
trate  farther  into  the  interstellar  medium,  enhancing  [C  ii]/ 
CO  (1-0)  and  [C  n]/[C  i]  (1-0)  and  increasing  the  beam  area 
filling  factor  of  [C  n]  relative  to  CO  and  [C  i]  (e.g.,  Maloney 
&  Black  1988;  Israel  et  al.  1996;  Madden  et  al.  1997; 
Madden  2000). 

The  peak  wavelength  of  the  dust  emission  from  high-U 
regions  will  depend  on  the  magnitude  of  U.  Such  regions 
can  help  explain  both  the  relatively  warm  60/100  /zm  colors 
in  ULIRGs  (e.g.,  Sanders  et  al.  1988)  and  the  weakness  of 
the  PAH  emission  relative  to  the  total  FIR  emission  as 
traced  by  the  6.2  /mi  PAH  feature  and  the  IRAS  12  /mi  flux. 
Boulanger  et  al.  (1988)  show  that  the  60  /my'  1 00  /mi  IRAS 
band  ratio  decreases  and  the  12  //m/25  //m  band  ratio 
increases  sharply  with  distance  from  the  central  star  within 
the  California  nebula  H  ii  region.  This  difference  in  behavior 
between  the  colors  within  H  ii  regions  and  in  the  surround¬ 
ing  PDR  is  virtually  identical  to  the  difference  in  both  the  60 
//m/100  //m  and  12  //m/25  /zm  colors  of  ULIRGs  compared 
with  less  luminous  normal  and  starburst  galaxies,  highly 
suggestive  of  an  increase  in  the  contribution  from  dust  in 
ionized  gas  to  the  FIR  emission  in  ULIRGs.  In  the 
California  nebula,  the  60  //m  flux  is  relatively  constant 
throughout  the  H  ii  region,  while  the  12  /zm  emission  has  a 
strong  peak  at  the  edge  of  the  H  n  region.  These  effects  are 
attributed  to  the  destruction  of  the  12  //m  emitters  and  the 
presence  of  large  grains  throughout,  because  the  emission 
from  large  grains  gradually  shifts  toward  shorter  wave¬ 
lengths  as  the  equilibrium  temperature  increases.  Mouri, 
Kawara,  &  Taniguchi  (1997)  argue  that  the  mid-lR  emis¬ 
sion  in  starburst  galaxies  does  in  fact  come  from  large  grains 
within  H  ii  regions.  Boulanger  et  al.  attribute  the  destruc¬ 
tion  of  the  12  /zm  emitters  to  the  increased  strength  of  the 
radiation  field  closer  to  the  star.  Desert,  Boulanger,  &  Puget 
(1990)  model  this  behavior  in  more  detail  with  a  three- 
component  model,  but  state  that  the  total  intensity  of  the 
radiation  field  is  a  less  critical  factor  than  the  presence  of 
hard  UV  photons.  Hard  UV  photons  are  present  within  H  ii 
regions  and,  of  course,  in  ionized  regions  around  AGNs.  On 
the  other  hand,  in  PDR  environments,  hard  UV  photons 
are  largely  absent,  and  for  the  G0  values  derived  here  for 
ULIRGs,  destruction  of  PAHs  is  not  likely  to  occur. 

Is  it  reasonable  that  high-U  conditions  contribute  in  an 
energetically  important  way  in  ULIRGs,  as  we  suggest? 
Using  mid-IR  spectroscopy  to  estimate  the  ratio  of  the  ion¬ 
izing  luminosity  ULyc  to  the  bolometric  luminosity  Ubob 
Genzel  et  al.  (1998)  find  ULyc/Ubol  values  in  ULIRGs  to  be 
nearly  half  that  in  their  starburst  templates.  Systematic 
beam  size  effects  in  their  starburst  templates  could  possibly 
make  the  discrepancy  larger  by  a  significant  factor.  In  our 
scenario,  measured  Q/L  values  would  be  significantly  lower 
in  ULIRGs  if  their  intrinsic  values  are  similar.  Given  the 
uncertainties  in  the  extinction  estimates  and  possible  beam 
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size  effects  in  the  nearby  starburst  nuclei,  further  measure¬ 
ments  are  warranted. 

Why  might  high -U  [=<b(H )/nec\  conditions  exist  in 
ULfRGs  and  what  possible  configurations  might  we 
envisage?  Given  the  evidence  for  relatively  higher  densities 
in  ULIRGs,  the  high-C/  conditions  likely  arise  from  high 
<1>(H)  due  to  high  Lyman  continuum  luminosities  impinging 
on  material  at  small  radii,  rather  than  from  low  electron 
densities.  These  conditions  could  arise  in  highly  compact 
H  n  regions,  or  they  could  arise  from  differences  in  the 
global  conditions  in  ULIRGs  compared  to  luminous  IR  gal¬ 
axies,  such  as  might  be  present  in  very  compact  luminous 
starburst  regions  where  the  ionized  regions  of  large  numbers 
of  individual  H  n  regions  merge  and  are  possibly  sur¬ 
rounded  by  a  global  PDR.  Here,  we  might  expect  high  U  to 
be  accompanied  by  high  G0. 

Lor  the  most  extreme  high- U  cases,  such  as  might  be 
present  in  buried  or  partially  obscured  AGNs,  the  dust 
emission  at  the  dust  sublimation  radius  must  peak  at  near- 
IR  wavelengths  (Dudley  &  Wynn-Williams  1997).  In  this 
case,  radiative  transfer  effects  could  produce  excess  LIR 
emission  arising  from  a  central  source  that  is  not  associated 
with  the  normal  PDR  tracers  in  these  galaxies  and  may  suf¬ 
fer  higher  extinction  than  the  PDR  components.  This 
explanation  may  be  responsible  for  the  observed  LIR 
molecular  absorption  seen  toward  ULIRGs  (Lischer  et  al. 
1999).  Such  entities  could  be  present  together  with  a  compo¬ 
nent  of  more  normal  and  possibly  less  obscured  PDR  and 
starburst  regions.  These  PDR  and  starburst  regions  could 
dominate  the  line  emission  and  thus  mask  the  presence  of  an 
obscured  central  source.  In  this  paper,  we  do  not  try  to 
ascertain  quantitatively  which  of  these  scenarios  is 
dominant,  although  pioneering  work  along  these  lines  by 
Laurent  et  al.  (2000)  has  attempted  to  quantify  separate 
PDR,  H  ii  region,  and  AGN  contributions  to  the  mid-IR 
spectral  signatures  of  galaxies.  Rather,  one  of  the  aims  of 
this  paper  is  to  quantify  the  observational  evidence  that 
high-  U  effects  are  important  and  to  emphasize  the  need  for 
further  observational  and  modeling  efforts. 

3.6.  Implications  for  More  Distant  Sources 

As  discussed  in  Paper  I,  an  important  implication  of  the 
[C  ii]  deficit  in  ULIRGs  is  that  the  usefulness  of  the  [C  ii]  line 
in  cosmological  work  could  be  diminished  if  cosmological 
sources  resemble  this  local  sample.  As  a  result  of  the  present 
work,  we  can  begin  to  provide  an  indication  of  what  this 
means  quantitatively.  In  the  present  sample,  only  3  of  15 
ULIRGs  resemble  normal  and  starburst  galaxies  in  the 
ratio  of  [C  n]-to-dust  continuum,  so  that  an  experiment  that 
could  just  detect  a  cosmological  ultraluminous  IR  source, 
such  as  those  detected  in  continuum  using  SCUBA  with 
normal  [C  ii]  properties,  would  need  to  observe  five  sources 
for  one  [C  ii]  detection.  Specifically,  for  the  high-redshift 
sources  considered  by  Stark  (1997),  the  per  season  detection 
rate  estimate  would  be  reduced  from  200  to  40.  Such  a  detec¬ 
tion  rate  should  still  provide  useful  information  for  under¬ 
standing  star  formation  as  a  function  of  redshift  (e.g., 
Madau,  Pozzetti,  &  Dickinson  1998),  since  lower  limits  to 
the  massive  star  formation  rate  can  be  calculated  directly 
from  the  measured  [C  ii]  intensity.  Given  the  presence  of 
upper  limits  to  the  [C  ii]  flux  in  our  sample,  it  is  not  possible 
to  say  what  would  be  required  to  detect  the  entire  sample  at 
cosmological  distances,  but  it  is  clear  that  should  the  [C  n] 


Fig.  9. — Ratio  of  [C  n]  intensity  to  FIR  flux  plotted  as  a  function  of  IR 
luminosity  for  the  present  sample  ( circles )  and  for  those  galaxies  in  the 
Malhotra  et  al.  (2001)  sample  ( diamonds )  that  appear  in  the  IRAS  Bright 
Galaxy  Sample.  Sources  plotted  as  open  symbols  may  have  both  con¬ 
tinuum  and  line  emission  that  extend  beyond  the  LWS  beam  (see  text).  The 
sample  medians  calculated  (ignoring  upper  limits)  for  the  filled  diamond 
symbol  galaxies  and  for  the  filled  circle  galaxies  with  Ljr  greater  than  or 
equal  to  1012  L0  are  plotted  by  solid  lines. 

deficit  hold,  it  may  become  a  driver  for  more  ambitious 
future  surveys. 

In  Ligure  9,  the  ratio  of  [C  ii]  to  LIR  flux  is  plotted  against 
8-1000  /in i  IR  luminosity  for  both  the  sample  considered 
here  and  the  Malhotra  et  al.  (2001)  sample.  Lor  galaxies  not 
listed  in  Table  2,  we  take  the  IR  luminosities  from  the  IRAS 
Bright  Galaxy  Sample  (Soifer  et  al.  1989)  or  its  extension 
(Sanders  et  al.  1995).  As  a  precaution,  we  have  plotted  as 
open  symbols  sources  that  are  listed  as  resolved  at  25  or  60 
fan  or  both,  although  we  have  made  use  of  the  LWS  data  to 
estimate  the  LIR  flux  where  appropriate  in  the  present  sam¬ 
ple.  It  is  of  interest  that  for  the  strongest  upper  limit  in  this 
plot  (NGC  4418),  the  power  source  of  the  LIR  is  considered 
to  be  a  buried  AGN  based  on  its  mid-IR  spectral  properties 
(Roche  et  al.  1986;  Dudley  &  Wynn-Williams  1997;  Spoon 
et  al.  2001).  The  data  in  Ligure  9  suggest  a  break  occurring 
near  1012  Le,  although  additional  data,  say  with  the 
Herschel  and  Astro-F  space  telescopes,  are  required  to  rule 
out  a  continuing  decline  at  higher  luminosities  and  to  better 
define  the  sharpness  of  the  break. 

Since  the  [C  n]  158  /mi  line  is  one  of  the  brightest  lines  in 
normal  and  starburst  galaxies  (and  ULIRGs)  in  the  wave¬ 
length  range  where  dust  extinction  plays  a  lesser  role,  its 
usefulness  as  an  astrophysical  diagnostic  remains  undimin¬ 
ished  notwithstanding  the  cautionary  implications  of  the 
present  work.  In  the  context  of  the  high-  U  hypothesis  dis¬ 
cussed  above,  sources  with  low  dust  abundance  relative  to 
hydrogen  would  have  to  achieve  an  even  higher  U  to  pro¬ 
duce  a  similar  diminution  of  the  [C  ii]  line  intensity  relative 
to  the  continuum.  Lor  such  scenarios,  we  might  expect  the 
break  suggested  in  Ligure  9  to  migrate  to  higher  luminosi¬ 
ties  if  U  and  L  are  generically  correlated  and  enrichment  is  a 
gradual  function  of  cosmic  time  in  individual  sources.  If 
such  is  the  case,  then  mid-  and  far-infrared  fine-structure 
lines,  which  are  preferred  ISM  abundance  diagnostics  (e.g., 
Simpson  et  al.  1998,  1995),  as  well  as  the  [C  ii]  158  //m  line 
itself,  will  be  easily  observed  in  moderate-luminosity  cosmo¬ 
logical  ultraluminous  IR  sources  and  can  provide  a  check 
on  this  scenario.  These  could  be  the  majority  of  such  sources 
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if  the  luminosity  function  slope  near  1012  L...  at  earlier  times 
is  similar  to  its  present  value  (see,  e.g.,  Soifer  &  Neugebauer 
1991;  Kim  &  Sanders  1998). 

Alternatively,  since  ULIRGs  are  mainly  associated  with 
major  mergers,  it  may  be  that  currently  known  cosmological 
ultraluminous  IR  sources  will  also  be  associated  with  (pairs 
of)  well-formed  and  thus  almost  certainly  pre-enriched 
systems.  If  this  is  the  case,  the  break  near  1012  Le  may  be 
intrinsic  to  the  properties  involved  in  the  formation  proc¬ 
esses  associated  with  ULIRGs.  In  this  scenario,  migration 
of  the  break  in  Figure  9  may  not  be  observable  at  submillim¬ 
eter  wavelengths.  A  more  detailed  discussion  of  the  cosmo¬ 
logical  implications  of  the  break  in  the  [C  n]  emission 
around  1012  LQ  in  light  of  the  above  scenarios  can  be  found 
in  Dudley  et  al.  (2003). 

4.  SUMMARY 

In  Paper  I,  we  reported  the  detection  of  a  deficiency  in  the 
[C  ii]  158  /mi  emission  relative  to  the  FIR  continuum  in  a 
small  sample  of  ULIRGs.  In  this  study,  we  combine  all  the 
available  [C  ii]  spectra  obtained  with  ISO  to  enlarge  the 
sample  from  7  to  15  ULIRGs,  12  of  which  show  detectable 
[C  ii]  emission.  We  find  that  the  enlarged  sample  of  ULIRGs 
confirms  our  earlier  observation  that  the  [C  ii]  flux  in 
ULIRGs  is  about  an  order  of  magnitude  lower  than  that 
seen  from  nearby  normal  and  starburst  galaxies.  As  dis¬ 
cussed  in  Paper  I,  this  result  could  have  important  implica¬ 
tions  for  the  use  of  the  [C  ii]  line  as  an  eventual  tracer  of 
high-redshift  protogalaxies,  which  may  resemble  the 
disturbed  ULIRGs  at  low  redshift  observed  in  this  study. 

In  an  attempt  to  address  the  origin  of  the  [C  ii]  deficit,  we 
combine  the  [C  ii]  data  with  ISO  LWS  spectra  of  the  [O  i]  63 
fim  line  fluxes  and  [O  i]  145  /an  upper  limits  in  the  brightest 
galaxies  in  our  sample;  like  [C  ii],  the  [O  i]  lines  are  thought 
to  arise  primarily  in  PDR  gas.  We  have  combined  the  [O  i] 
measurements  with  other  recent  observations  of  the  atomic 
fine-structure  and  molecular  rotation  line  radiation  and 
PAH  emission  in  ULIRGs  to  explore  the  origin  of  the  [C  ii] 
deficit.  In  particular,  we  have  for  the  first  time  examined  in 
detail  the  self-consistency  and  applicability  of  the  PDR 
assumptions  in  the  context  of  a  filling  factor  analysis  for 
many  of  the  most  often  observed  diagnostics. 

Based  on  the  [C  ii]  and  [O  i]  observations  alone,  we  find 
that  PDR  models  suggest  that  ULIRGs  are  characterized 
by  high  Go/n  PDRs,  but  only  if  the  PDR  model  assumption 
that  the  [C  ii],  [O  i],  and  FIR  emission  components  are 
coextensive  holds,  i.e.,  both  the  lines  and  the  FIR  trace  PDR 
gas  and  have  identical  beam  area  filling  factors.  Within  the 
context  of  these  models  the  observed  decrease  in  [C  ii]/FIR 
and  6.2  pm  PAH/FIR  with  60  /im/100  pm  (which  scales 
with  G0  in  the  PDR)  provides  apparent  support  for  the  pres¬ 
ence  of  high  Go/n  PDRs  in  ULIRGs.  We  do  find,  however, 
that  any  PDR  model  that  can  account  for  the  observed 
[C  ii],  [O  i]  63  and  145  pm,  and  FIR  emission  does  not 
compare  well  with  previous  observations  of  the  neutral/ 
molecular  ISM  in  ULIRGs. 

We  use  PDR  models  to  derive  beam  area  filling  factors 
for  [C  ii],  CO  (1-0),  and  [C  i]  (1-0).  Inconsistencies  in  the 
derived  beam  area  filling  factors  suggest  that  the  assump¬ 
tion  that  the  [C  ii]  and  FIR  fluxes  arise  from  the  same  region 
does  not  apply  for  ULIRGs.  If  true,  we  argue  that  non- 
PDR  components  that  dominate  the  FIR  emission  but  do 
not  contribute  to  the  [C  n]  flux  can  account  for  the  [C  n] 


deficit.  (Beam  filling  factor  effects  do  not  rule  out  an  added 
contribution  from  high  Go/n  PDRs  to  the  [C  n]  deficit;  they 
only  preclude  the  use  of  one-dimensional,  plane-parallel 
models  as  a  means  to  justify  a  high  Go/n  scenario.)  We  point 
out  that,  for  previous  PDR  analyses  that  assume  that  the 
FIR  flux  arises  solely  from  PDRs,  disparate  PDR  physical 
conditions  are  derived  depending  on  the  lines  used:  [C  n] 
and  [O  i]  alone  yield  high  G0/n  solutions,  whereas  when  CO 
(1-0)  and  [C  i]  (1-0)  are  compared  with  [C  n],  high  n  and 
high  G0  are  derived.  This  discrepancy  further  suggests  that 
non-PDR  contributors  to  the  FIR  emission  are  at  work — a 
finding  supported  by  the  observed  similarity  in  the  [C  ii]  to 
6.2  /im  PAH  ratio  in  ULIRGs  and  less  luminous  galaxies. 
In  addition,  we  conclude  that  soft  UV  radiation  fields  aris¬ 
ing  from  aging  starbursts,  optical  depth  and  FIR  extinction 
effects,  and  high-/;  and  high-G0  (n  >  105  cm~3,  Go  >  104) 
PDRs  are  not  compatible  with  the  full  set  of  FIR  line  and 
continuum  observations  (although  mid-IR  extinction  may 
play  a  role  in  the  PAH  deficit  relative  to  FIR). 

As  an  example  of  a  non-PDR  FIR  component,  enhanced 
dust  absorption  in  the  ionized  gas  may  be  able  to  account 
for  the  implied  beam  area  filling  factor  differences  between 
[C  ii]  and  CO,  [C  i],  and  FIR.  Enhanced  dust  absorption  of 
ionizing  photons  in  ionized  regions  could  give  rise  to  signifi¬ 
cant  FIR  emission  from  the  photoionized  gas.  This  same 
gas  would  not  produce  strong  [C  n]  emission.  Such  an  effect 
could  give  rise  to  smaller  beam  area  filling  factors  for  [C  ii] 
than  for  FIR  and  produce  a  [C  ii]  deficit  relative  to  FIR  as 
observed.  Photoionized  gas  characterized  by  a  high  U  can 
produce  such  enhanced  dust  absorption  and  excess  FIR 
emission.  High-  U  ionized  regions  also  provide  an  explana¬ 
tion  that  simultaneously  accounts  for  the  suppressed  FIR 
fine-structure  line  and  PAH  emission  and  the  warmer  FIR 
colors  and  dust  temperatures  of  ULIRGs  compared  to  less 
luminous  galaxies. 

Should  the  [C  ii]  deficit  hold  in  high-redshift  protogalax¬ 
ies,  it  could  have  an  impact  on  the  use  of  the  [C  ii]  line  for 
future  surveys.  An  estimate  of  the  detection  rate  of  such 
objects  based  on  the  [C  ii]  observations  presented  here 
would  imply,  however,  that  future  [C  n]  observations  may 
still  yield  useful  information  for  understanding  star  forma¬ 
tion  as  a  function  of  redshift.  Here,  the  dust  abundance  of 
these  high-redshift  objects  will  play  an  important  role,  since 
sources  with  low  dust  abundance  relative  to  hydrogen 
would  have  to  achieve  an  even  higher  U  to  produce  a  [C  ii] 
deficit  similar  to  that  seen  for  the  ULIRGs.  Future  observa¬ 
tions  of  the  [C  ii]  line  in  high-redshift  sources  may  therefore 
help  to  provide  a  check  on  the  proposed  high-  U  scenario  in 
ULIRGs. 
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APPENDIX 

In  this  appendix,  we  discuss  the  basis  for  the  assumption  in  this  work  that  the  [O  i]  and  [C  n]  lines  in  extragalactic  sources 
are  optically  thick  and  thus  arise  from  the  near  side  of  the  clouds,  as  adopted  in  the  Kaufman  et  al.  (1999)  PDR  models.  For 
plane-parallel  models,  this  produces  a  factor  of  0.5  in  the  ratios  of  these  line  fluxes  to  the  total  FIR  flux,  as  discussed  in  §  3. 1 . 

Kaufman  et  al.  (1999)  model  a  typical  galactic  PDR  associated  with  a  star-forming  region  in  which  the  line  width  is  1.5  km 
s_1  and  Ay^lO,  corresponding  to  a  total  column^  x  1022  cm~2.  This  analysis  can  be  applied  to  galaxies  as  long  as  most  PDRs 
along  the  same  line  of  sight  have  non-overlapping  projected  velocities.  Of  the  total  column,  the  [O  i]  emitting  region  on  the  far 
side  of  a  cloud  corresponds  to  approximately  Ay  =  3,  or  a  column  of  N  =  77(14  i)  +  2A(H2)  =  6  x  1021  cm~2.  Thus,  the  [O  i] 
line  emission  passes  through  a  column  of  N  =  1.4  x  1022  (from  the  far  side  of  the  cloud  to  the  observer).  In  the  cloud  center,  of 
extent  approximately  Ay  =  4,  the  O  is  partly  tied  up  in  CO,  with  the  remainder  mainly  in  the  form  of  atomic  O.  From  Table  4 
of  Tielens  &  Hollenbach  (1985),  an  optical  depth  of  1  occurs  in  the  [O  i]  line  for  a  column  of  only  N  =  7.8  x  1020  ciru2.  This 
column  should  be  revised  for  the  differences  in  gas-phase  abundance  (5  x  10~4  and  3  x  10  4)  and  line  width  (2.7  and  1.5  km 
s  1 )  between  that  of  Tielens  &  Hollenbach  and  Kaufman  et  al.  (1999),  respectively,  using  only  the  atomic  O  which  is  not  tied 
up  in  CO  (1 .6  x  10  4)  in  the  central  region.  Taking  these  differences  into  account,  estimates  of  the  optical  depth  in  the  [O  i]  line 
corresponding  to  these  near  side  and  central  parts  of  the  column  are: 

near  side  (Ay  =  3), 

6  x  1021  cm-2  3  x  10~4  2.7  km  s^1  Q 

^63 /an  “  ?  g  x  1Q20  cm-2  5  x  1Q-4  1.5  km  S"1  =  ^  ’ 

central  region  (Av  =  4), 

8  x  1021  cm~2  1.6  x  10-4  2.7  km  s_1 

*Wm  =  7  g  x  1020cm-2  5x10-4  i.5kms-1=  '  ' 

This  yields  a  total  optical  depth  of  r[0  i]63  =  14.2.  With  an  optical  depth  for  [O  i]  63  /tm  on  the  order  of  14,  a  fraction  of  only 

e~  14  «  10  6  of  the  emission  penetrates  from  the  far  side  through  the  near  side  of  the  cloud.  (We  note  that  the  fact  that  self¬ 
absorption  is  rarely  seen  means  that  there  is  little  cold  gas  between  the  last  emitting  region,  the  near  side  of  the  cloud,  and  us, 
not  that  the  line  is  optically  thin.  In  addition,  the  temperature  structure  of  the  PDR,  with  hotter  regions  near  the  surface  and 
cooler  regions  in  the  interior,  is  not  expected  to  produce  self-absorbed  profiles  even  for  optically  thick  lines  emitted  from  the 
near  side  of  the  cloud.) 

For  C+,  Tielens  &  Hollenbach  (1985)  calculate  an  optical  depth  of  1  at  a  column  density  of  1.2  x  1021  cm~2.  This  column 
density  should  be  revised  by  a  factor  of  (3.0  x  10  4/l .4  x  10~4)(1.5  km  s_1/2.7  km  s-1)  =  1 .2  to  yield  a  column  of  1 .4  x  1021 
ciru2  required  to  produce  an  optical  depth  of  1 .  The  C+  column  extends  for  an  A  v  of  about  2  into  the  PDR,  but  the  optical 
depth  is  dependent  on  G{)  and  n  (and  T).  An  A  [/  of  2  would  give  an  optical  depth  in  the  line  of  about  2.8.  From  our  detailed 
PDR  models,  at  the  typical  values  Go  =  103  and  n  =  103,  we  find  an  optical  depth  of  about  1  in  [C  ii]  158  /an.  Thus,  we  expect 
that  less  than  40%  of  the  [C  ii]  emission  can  penetrate  from  the  far  side  of  the  cloud  through  the  near  side  to  us,  i.e.,  the  [C  ii] 
emission  is  marginally  thick.  We  note  that  the  fact  that  [C  ii]  absorption  is  seen  toward  the  Galactic  center  at  some  velocities 
(Vastel  et  al.  2002)  arises  from  a  unique  perspective  in  which  we  view  many  clouds  that  are  predominantly  moving  tangentially 
to  our  line  of  sight.  This  produces  a  coherent  path  length  of  [C  ii]  emission,  which  leads  to  self-absorption.  The  same 
perspective  does  not  generally  apply  when  we  view  an  external  galaxy. 
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